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We review recent progress of modified gravity models of dark energy-based on f(R) gravity, 
scalar-tensor theories, braneworld gravity, Galileon gravity, and other theories. In f(R) gravity 
it is possible to design viable models consistent with local gravity constraints under a chameleon 
mechanism, while satisfying conditions for the cosmological viability. We also construct a class 
of scalar-tensor dark energy models based on Brans-Dicke theory in the presence of a scalar-field 
potential with a large coupling strength Q between the field and non-relativistic matter in the 
Einstein frame. We study the evolution of matter density perturbations in f(R) and Brans-Dicke 
theories to place observational constraints on model parameters from the power spectra of galaxy 
clustering and Cosmic Microwave Background (CMB). 

The Dvali-Gabadazde-Porrati braneworld model can be compatible with local gravity constraints 
through a nonlinear field self-interaction □0(9 M 09 M 0) arising from a brane-bending mode, but the 
self-accelerating solution contains a ghost mode in addition to the tension with the combined data 
analysis of Supernovae la (SN la) and Baryon Acoustic Oscillations (BAO). The extension of the 
field self-interaction to more general forms satisfying the Galilean symmetry d^cj) — > d^cj) + 6 M in 
the flat space-time allows a possibility to avoid the appearance of ghosts and instabilities, while the 
late-time cosmic acceleration can be realized by the field kinetic energy. We study observational 
constraints on such Galileon models by using the data of SN la, BAO, and CMB shift parameters. 

We also briefly review other modified gravitational models of dark energy-such as those based on 
Gauss-Bonnet gravity and Lorentz- violating theories. 
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I. INTRODUCTION 



The cosmic acceleration today has been supported by independent observational data such as the Supernovae-type 
la (SN la) [H, Hj], the Cosmic Microwave Background (CMB) temperature anisotropies measured by WMAP plli}, 
and Baryon Acoustic Oscillations [U @ . The origin of dark energy responsible for this cosmic acceleration is one of 
the most serious problems in modern cosmology |7H2fj|. The cosmological constant is one of the simplest candidates 
for dark energy, but it is plagued by a severe energy scale problem if it originates from the vacuum energy appearing 
in particle physics [2l|. 

The first step toward understanding the nature of dark energy is to clarify whether it is a simple cosmological 
constant or it originates from other sources that dynamically change in time. The dynamical models can be distin- 
guished from the cosmological constant by considering the evolution of the equation of state of dark energy (= u>de)- 
The scalar field models of dark energy such as quintessence [22ll26| and k-essence [13, [H| predict a wide variety of 
the variation of wde> but still the current observational data are not sufficient to distinguish those models from the 
A-Cold-Dark-Matter (ACDM) model. Moreover it is generally difficult to construct viable scalar-field models in the 
framework of particle physics because of a very tiny mass (m^ < 10 -33 eV) required for the cosmic acceleration today 
[illlj. 

There exists another class of dynamical dark energy models based on the large-distance modification of gravity. The 
models that belong to this class are f(R) gravity J3l|-[3J| (/ is function of the Ricci scalar R), scalar-tensor theories 
[35M40l | , braneworld models [4l[ , Galileon gravity |42| . Gauss-Bonnet gravity [43|, , and so on. An attractive feature 
of these models is that the cosmic acceleration can be realized without recourse to a dark energy matter component. 
If we modify gravity from General Relativity (GR), however, there are tight constraints coming from local gravity 
tests as well as a number of observational constraints. Hence the restriction on modified gravity models is in general 
stringent compared to modified matter models (such as quintessence and k-essence). 

For example, an f (R) model of the form f(R) = R — fi 2 ( n+1 ^ / R n (n > 0) was proposed to explain the late-time 
cosmic acceleration [§2|, 33:] (see also Refs. |45l - [5Cl | for early works). However this model suffers from a number of 
problems such as the incompatibility with local gravity constraints [5ll - [54j |, the instability of density perturbations 
[55M58I ] , and the absence of a matter-dominated epoch [5^, [6(| • As we will see in this review there are a number of 
conditions required for the viability of f(R) dark energy models [55l - [58ll6lT[67| . which stimulated to propose viable 
models [6814721 ]. 

The simplest version of scalar-tensor theories is so-called Brans-Dicke theory in which a scalar field <p couples to 
the Ricci scalar R with the Lagrangian density C — ipR/2 — (wbd/2</?) (Vc/?) 2 , where wbd is a so-called Brans-Dicke 
parameter [73J . GR can be recovered by taking the limit wbd — >• 00. If we allow the presence of the field potential 
U(tp) in Brans-Dicke theory, f(R) theory in the metric formalism is equivalent to this generalized Brans-Dicke theory 
with the parameter wbd = [5l|, [7J| . By transforming the action in generalized Brans-Dicke theory ( "Jordan frame" ) 
to an "Einstein frame" action by a conformal transformation, the theory in the Einstein frame is equivalent to a 
coupled quintessence scenario [75|] with a constant coupling Q satisfying the relation l/(2Q 2 ) = 3 + 2wbd UM- For 
example, f(R) theory in the metric formalism corresponds to the constant coupling Q = — l/y6, i.e. wbd — 0. For 
\Q\ of the order of unity it is generally difficult to satisfy local gravity constraints unless some mechanism can be at 
work to suppress the propagation of the fifth force between the field and non-relativistic matter. It is possible for 
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such large-coupling models to be consistent with local gravity constraints [U [6|| [77l - [80| through the so-called 
chameleon mechanism |8lL l82j. provided that a spherically symmetric body has a thin-shell around its surface. 

A braneworld model of dark energy was proposed by Dvali, Gabadadze, and Porrati (DGP) by embedding a 3- 
brane in the 5-dimensional Minkowski bulk spacetime [4l|. In this scenario the gravitational leakage to the extra 
dimension leads to a self-acceleration of the Universe on the 3-brane. Moreover a longitudinal graviton (i.e. a brane- 
bending mode (j>) gives rise to a nonlinear self-interaction of the form (r^/m p \) □0(9 M (/><9 AI </>) through the mixing with a 
transverse graviton, where r c is a cross-over scale (of the order of the Hubble radius Hq 1 today) and m v \ is the Planck 
mass [83l. l84j. In the local region where the energy density p is much larger than r~ 2 m£\ the nonlinear self-interaction 
can lead to the decoupling of the field from matter through the so-called Vainshtein mechanism [85], which allows a 
possibility for the consistency with local gravity constraints. However the DGP model suffers from a ghost problem 
|86l - l88j , in addition to the difficulty for satisfying the combined observational constraints of SN la and B AO [89|-[94[ ■ 

The equations of motion following from the self-interacting Lagrangian □0(<9''</><9 A1 </>) present in the DGP model are 
invariant under the Galilean shift d^tfi — > d^ + b^ in the Minkowski background. While the DGP model is plagued by 
the ghost problem, the extension of the field self-interaction to more general forms satisfying the Galilean symmetry 
may allow us to avoid the appearance of ghosts. Nicolis et al. [42| showed that there are only five field Lagrangians 
d (i = 1, • • • , 5) that respect the Galilean symmetry in the Minkowski background. In Refs. [§5|,[96| these terms were 
extended to covariant forms in the curved space-time. In addition one can keep the equations of motion up to the 
second-order, while recovering the Galileon Lagrangian in the limit of the Minkowski space-time. This property is 
welcome to avoid the appearance of an extra degree of freedom associated with ghosts. In fact, Refs. [ItIIsH derived 
the viable model parameter space in which the appearance of ghosts and instabilities associated with scalar and tensor 
perturbations can be avoided. Moreover the late-time cosmic acceleration is realized by the existence of a stable de 
Sitter solution. We shall review the cosmological dynamics of Galileon gravity as well as conditions for the avoidance 
of ghosts and instabilities. 

In order to distinguish between different models of dark energy based on modified gravitational theories, it is 
important to study the evolution of cosmological perturbations as well as the background expansion history of the 
Universe. In particular, the modified growth of matter perturbations S m relative to the ACDM mod el changes the 
matter power spectrum of large-scale structures (LSS) as well as the weak lensing spectrum [99|- jl26j . Moreover the 
modification of gravity manifests itself for the evolution of the effective gravitational potential ij> related with the 
Integrated-Sachs- Wolfe (ISW) effect in CMB anisotropies. We shall review a number of observational signatures for 
the modified gravitational models of dark energy. 

This review is organized as follows. In Sec.|TT]we construct viable dark energy models based on f(R) theories after 
discussing conditions for the cosmological viability as well as for the consistency with local gravity tests. In Sec. lIIII wc 
show that, in Brans-Dicke theories with large matter couplings, it is possible to design the field potential consistent 
with both cosmological and local gravity constraints. In Sec. IIVI we derive the field equations in the DGP model 
and confront the model with observations at the background level. In Sec. |V]we review the cosmological dynamics 
based on Galileon gravity as well as conditions for the avoidance of ghosts and Laplacian instabilities. In Sec. I VII we 
briefly mention other modified gravity models of dark energy based on Gauss-Bonnet gravity and Lorentz- violating 
theories. In Sec. IVIII we study observational signatures of dark energy models based on f(R) gravity, Brans-Dicke 
theory, DGP model, and Galileon gravity, in order to confront them with the observations of LSS, CMB, and weak 
lensing. Sec. I Villi is devoted to conclusions. 

Throughout the review we use the units such that c = H = hs = l> where c is the speed of light, h is reduced 
Planck's constant, and fee is Boltzmann's constant. We also adopt the metric signature (— , +, +, +). 



II. f(R) GRAVITY 



We start with the action in f(R) gravity: 



where k 2 = 8irG (G is a bare gravitational constant), g is a determinant of the metric g^u, f(R) is an arbitrary 
function in terms of the Ricci scalar R, and Cm is a matter action with matter fields "I'm- Variation of the action ([!]) 
with respect to g^ leads to the following field equation 

F{R)R llv {g) - ±f(R)g MV - V„V„F(fl) + g^UF{R) = k 2 T^ , (2) 
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where F(R) = fa = df /dR, R^ v is a Ricci tensor, and T^ v — — (2/ y/—g)SCM l&g^ v is an energy- momentum tensor 
of matter. The trace of Eq. ([2]) gives 

3ElF{R) + F{R)R-2f(R) = k 2 T, (3) 

where T — g^T^ = —p + 3P. Here p and P are the energy density and the pressure of matter, re spec tively. 

Regarding the variation of the action ([TJ , there is another approach called the Palatini formalism [l27j | in which g^ 
and the affine connection are treated as independent variables. 1 The resulting field equations are second-order 
[l32l - ll37l ] and the cosmological dynamics of dark energy models have been studied by a number of authors [l38l - ll45j | . 
However f(R) theory in the Palatini for malism gives rise to a large coupling between a scalar field degree of freedom 
and ordinary matter [135H1371 Il46l - |l50l | , which implies difficulty for compatibility with standard models of particle 
physics. This large coupling als o leads t o significant growth of matter density perturbations, unless the models are 
very close to the ACDM model [HMH. 

In the following we focus on the variational approach (so called the metric formalism) given above. The Einstein 
gravity without a cosmological constant corresponds to f{R) = R and F{R) — 1, so that the term DF(R) in Eq. ([3]) 
vanishes. Since in this case R = —k 2 T = K 2 (p — 3P), the Ricci scalar R is directly determined by matter. In f(R) 
gravity with a non-linear term in R, DF(R) does not vanish in Eq. <j3j) . Hence there is a propagating scalar degree of 
freedom, if) = F(R), dubbed "scalaron" in Ref. |l56l |. The trace equation ^ allows the dynamics of the scalar field 

i>. 

The de Sitter point corresponds to a vacuum solution with constant R. Since DF(R) — at this point, we obtain 

F(R)R-2f(R)=0. (4) 

Since the quadratic model f(R) — uR 2 satisfies this condition, it gives rise to an exact de Sitter solution. In the 
inflation model f(R) = R + aR 2 proposed by Starobinsky [l56l |. the accelerated cosmic expansion ends when the term 
aR 2 becomes s maller tha n the linear term R. It is possible to construct such f(R) inflation models in the framework 
of supergravity jl57l Il58l | . 

A. Cosmological dynamics in f(R) gravity 

We first study cosmological dynamics for the models based on f(R) theories in the metric formalism. In order to 
derive conditions for the cosmological viability of f(R) models we shall carry out general analysis without specifying 
the form of f(R). We consider a flat Friedmann-Lemaitre- Robertson- Walker (FLRW) background with the line 
element 

ds 2 = -dt 2 + a(i) 2 da; 2 , (5) 

where a(t) is a scale factor. For the matter Lagrangian Cm in Eq. (JT]) we take into account non-relativistic matter and 
radiation, whose energy densities p m and p r satisfy the usual continuity equations p m + 3Hp m — and p r + AH p r = 
respectively. Here H = a/a is the Hubble parameter and a dot represents a derivative with respect to cosmic time t. 
From Eqs. ([2]) and ((3|) we obtain 

3FH 2 = K 2 (p m + p r ) + (FR- f)/2-3HF, (6) 
2FH = -k 2 [p rn + (4/3)p r ] — F + HF , (7) 

where the Ricci scalar is given by 

R = 6(2H 2 + H) . (8) 
Let us introduce the following dimensionless variables: 

Xl ~ HF' X2 ~ 6FH 2 ' X3 ~ m 2 ' Xi - 3FH 2 ' 1 ' 



1 We also note that there is another approach for the variational principle— known as the metric-affine formalism-in which the matter 
Lagrangian Cm depends not only on the metric but also on the connection Tg^ [l28Hl3l]l . 
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together with the density parameters 



2 

3FH 2 



1 — Xl — X 2 — X 3 — X4, 



£ 4 , 



DE 



Xi + X 2 + X3 . 



It is straightforward to derive the following equations [62] : 



dxi 
dN 
dx 2 
dN 
dx 3 
dN 

dX4 

dN 



— 1 — X 3 — 3X2 + X x — X\X 3 + X4 
X\X 3 



m 

xix 3 



x 2 (2x 3 - 4- xi) 
- 2x 3 (x 3 - 2) , 



where N = In a and 



m 



m 

-2x 3 X4 + X1X4 , 

dlnF Rf RR 



dlnP 
din/ 
dln-R 



f,R 
Rf,R 
f 



El 

X2 



(10) 

(11) 

(12) 
(13) 
(14) 

(15) 
(16) 



From Eq. (|16[) one can express R as a function of x 3 /x 2 - Since m is a function of R, it follows that m is a function 
of r, i.e., m = m(r). The ACDM model, /(i?) = i? — 2A, corresponds to to = 0. Then the quantity to characterizes 
the deviation from the ACDM model. 

The effective equation of state of the system is given by 



-1 



2H 



1 



(2x 3 - 1) 



3H 2 3' 

In the absence of radiation [xn — 0) the fixed points for the dynamical system (ITT|) - (IT41 are 



Pi 

P-2 
P 3 
P4 

P 5 

Ph 



(xi,x 2 ,x 3 ) 
(xi,x 2 ,x 3 ) 
(xi,x 2 ,x 3 ) 
(xi,x 2 ,x 3 ) 

(xi,x 2 ,x 3 ) 

(X 1 ,X 2 ,X 3 ) 



(0,-1,2), O m = 0, w ee = -l, 
(-1,0,0), n m = 2 1 w cff = l/3, 
(1,0,0), a„ = 0, uyr = l/3, 
(-4,5,0), O m = 0, w ee = l/3, 
3to 1 + 4to 1 + 4m 



l + m' 2(1 + m) 2 ' 2(1 + to) 
2(1 -to) 1-4to (1-4to)(1 + to) 



1 



1 + 2to ' to(1 + 2m)' 



to(1 + 2to) 



m(7 + 10to) 
2(1 + to) 2 

VL m = 0, w cS 



WeS 



l + m 
2 — 5m — 6to 2 
3m(l + 2m) 



(17) 

(18) 
(19) 
(20) 
(21) 

(22) 
(23) 



The points P5 and Pq are on the line m(r) = — r — 1 in the (r, m) plane. 

Only the point P5 can be responsible for the matter-dominated epoch (fl m ~ 1 and w c s — 0). This is realized 
provided m is close to 0. In the (r, to) plane the matter point P5 exists around (r, to) = (—1,0). Either the point Pi 
or Pq can lead to the late-time cosmic acceleration. The former corresponds to a de Sitter point (w c g — — 1) with 
r = —2, in which case the condition (j4]) is satisfied. Depending on the values of to, the point Pq can be responsible 
for the cosmic acceleration 62] . In the following we shall focus on the case in which the matter point P5 is followed 
by the de Sitter point Pi. 

The stability of the fixed points is known by considering small perturbations Sxi (i = 1, 2, 3) around them [62]. For 
the point P5 the eigenvalues for the 3x3 Jacobian matrix of perturbations are 



3(l + m' 5 ), 



-3to 5 ± y/m 5 (256m 3 5 + 160m 2 - 31m 5 - 16) 

4TO5(TO5 + 1) 



(24) 



where 7715 = m{r§) and m' 5 = ^(^5) with r& —1. In the limit 1 777,5 1 "C 1 the latter two eigenvalues reduce to 
—3/4 ± -y/— 1/777,5. The f(R) models with 777,5 < show a divergence of the eigenvalues as TO5 — > —0, in which case 
the system cannot remain for a long time around the point P5. For example the model f(R) = R — a/R n with 77 > 








Figure 1: Four trajectories in the (r,m) plane. Each trajectory corresponds to the models: (i) ACDM, (ii) f(R) = (R b — A) c , 
(iii) f(R) = R- aR n with a > 0, < n < 1, and (iv) m(r) = -C(r + l)(r 2 + ar + b). Here P M , Pa and P B are the matter 
point P$, the de Sitter point Pi, and the accelerated point P§, respectively. From Ref. |64l |. 



and a > falls into this category. On the other hand, if < m 5 < 0.327, the latter two eigenvalues in Eq. (l24l) arc 
complex with negative real parts. Then, provided that m' 5 > —1, the point P5 corresponds to a saddle point with 
a damped oscillation. Hence the Universe can evolve toward the point P5 from the radiation era and leave for the 
late-time acceleration. Then the condition for the existence of the saddle matter era is 

, . dm 

m(r)«+0, — > -l at r = -1 . 25) 
dr 

The first condition implies that the f(R) models need to be close to the ACDM model during the matter era. 
The eigenvalues for the Jacobian matrix of perturbations about the point Pi are 



3. _»± Vg^ , (26 , 



where mi = m(r = —2). This shows that the condition for the stability of the de Sitter point Pi is (62l . Il59l - ll61 | 



0<m(r=-2)<l. (27) 



The trajectories that start from the saddle matter point P5 with the condition (|25j) and then approach the stable de 
Sitter point Pi with the condition (|2~7| are cosmologically viable. 

Let us consider a couple of viable f(R) models in the (r, m) plane. The ACDM model, f(R) = R—2A, corresponds to 
m = 0, in which case the trajectory is a straight line from P5: (r, m) = (—1, 0) to Pi: (r, m) = (—2, 0). The trajectory 

(ii) in Fig.[T]represents the model f(R) = (R b —A) c [Hj], which corresponds to the straight line m(r) = [(1— c)/c]r+b— 1 
in the (r, m) plane. The existence of a saddle matter epoch requires the condition c > 1 and be ~ 1. The trajectory 

(iii) represents the model lolil ] 

f(R) = R-aR n (a > 0, < n < 1), (28) 

which corresponds to the curve m = nil + r)/r. The trajectory (iv) in Fig. Q] shows the model m(r) = — C(r + 1)(^ 2 + 
ar + 6), in which case the late-time accelerated attractor is the point Pq with (V3 — l)/2 < m < 1. 

In Ref. [62j it was shown that the variable m needs to be close to during the radiation-dominated epoch as well. 
Hence the viable f(R) models are close to the ACDM model, f(R) = R — 2A, in the region R 3> Ro (where Rq is 
the present cosmological Ricci scalar) . The Ricci scalar R given in Eq. ([5]) remains positive from the radiation era to 
the present epoch, as long as the it does not oscillate. As we will see in Sec. Ill Bl we require the condition /j? > 
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to avoid ghosts. Then the condition m > for the presence of the matter-dominated epoch translates to J^rr > 0. 
The model f(R) = R — a/R n (a > 0, n > 0) is not viable because the condition f.RR > is violated. We also note 
that the power-la w mo dels with f(R) — R n do not give rise to a successful cosmological trajectory [5^, [6(| (unlike 
the claims in Ref. |162| ). 

In order to derive the equation of state of dark energy to confront with SN la observations for the cosmologically 
viable models, we rewrite Eqs. ((6]) and as follows 2 : 

3AH 2 = k 2 ( Pm + Pr + p DE ) , (29) 
-2AH = n 2 [ Pm + (4/3)p r + pde + Pde] , (30) 

where A is some constant and 

k 2 Pde = (1/2)(FR- f)-3HF + 3H 2 (A-F), (31) 
k 2 P be = F + 2HF - (l/2)(FR- f) - (3H 2 + 2H)(A- F) . (32) 

Defining prjE and Pde in this way, one can show that these satisfy the usual continuity equation 

PDE + 3 J H'(p DE + J FbE) = 0. (33) 

The dark energy equation of state related with SN la observations is given by w^e = Pde/ Pde- From Eqs. (|2"9")l 
and (|30|) it follows that 

= 2AH + 3AH 2 + n 2 p r /3 w cS 

WDE 3 AH 2 - K 2 (p m + Pr ) - 1 - (F/A)n m ' [ ' 

where the last approximate equality in Eq. (|34j) is valid in the regime where the radiation density p r is negligible 
relative to the matter density. The viable f(R) models approach the ACDM model in the past, i.e. F — > 1 as R — » oo. 
In order to reproduce the standard matter era for the redshifts z 3> 1, one can choose A = 1 in Eqs. P9^) and (|30[) . 
Another possible choice is A — Fq, where Fq is the present value of F. Th is choice is suitable if the deviation of Fq 
from 1 is small (as in the scalar-tensor theory with a massless scalar field 164 - 166 ). In both cases the equation of 



167H169j . Thus f(R) gravity models 



164 

state wde can be smaller than —1 before reaching the de Sitter attractor [64 .171. 
give rise to a phantom equation of state without violating stability conditions of the system. 

As we see in Eq. (|34l) . the presence of non-relativistic matter is important to lead to the apparent phantom behavior. 
We wish to stress here that for viable f(R) models constructed to satisfy all required conditions [such as the models 
(|4"5j) , (|4*6"|) , and (14*51) we will discuss later] the ghosts are not present even if w de < — 1- A number of authors proposed 
some models to realize wde < — 1 without including non-relativistic matter [l70l Il7l| , which means that u>de = w e e 
from Eq. (|34| . However, such f(R) models usually imply the presence of ghosts 3 , because w e g < —1 corresponds to 
H > 0. 

The observational constraints on specific f(R) models ha ve been carried out in Refs. |172| - |175| from the background 
expansion history of the Universe (see also Refs. |176l 4l80] for the reconstruction of f(R) mo dels f rom observations). 
Since the deviation of wde from that in the ACDM model (i«de = — 1) is not so significant 68, 167], the viable models 
such as (|4"5j) . (|4"6")l . and (|48p can be consistent with the data fairly easily. In other words we do not obtain very tight 
bounds on model parameters from the information of the background expansion history only. However, the models 
can be more strongly constrained at the level of perturbations, as we will see in Sec. IVII Al 



B. Conditions for the avoidance of ghosts and tachyonic instabilities 

In this subsection we shall derive conditions for the avoidance of ghosts and tachyonic instabilities in f(R) theories. 
In doing so we expand the action ([1} up to the second-order by considering the following perturbed metric about the 
FLRW background 

ds 2 = -(1 + 2a)dt 2 - 2a(t)d^dtdx l + a 2 {t) {5 l] + 2^S %3 + 2d l d jl ) dx l dx 3 , (35) 



2 If the field equations are written in this form, we can also show that the background cosmological dynamics has a correspondence with 
equilibrium thermodynamics on the apparent horizon Il63(| . 

3 If a late-time de Sitter solution is a stable spiral, it happens that iode oscillates around —1 with a small amplitude, even for viable f(R) 
models. Here we are discussing the real ghosts out of this regime. 
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where a, /?, ip, 7 are scalar metric perturbations 181] . 

Introducing the perturbation SF for the quantity F — df/dR, one can construct the gauge- invariant curvature 
perturbation 



n = tb--^SF. 
F 



(36) 



Expandi ng the action ([1} without the matter source, we obtain the second-order action for the curvature perturbation 



<^ 2 



dt d 3 xa 3 Q s 



(37) 



where 



3F 2 



2k 2 F[H + F/(2F)} 2 



(38) 



The negative sign of Q s corresponds to a ghost field because of the negative kinetic energy. Hence the condition for 
the avoidance of ghosts is given by 



F > 0. 



(39) 



For the matte r sec tor the ghost does not appear for pm(l + wm)/wm > (where wm is the equation of state for the 
matter fluid) 183], which is satisfied for radiation (% = 1/3) and non-relativistic matter (wm — +0). 

If Q s is positive, the action (|37l) can be written in the following form by introducing the new variables u = z s lZ 
and z s = ay/Q s : 



8SW 



drj d 3 x 



1 I z" 



(40) 



where a prime represents a derivative with respect to the conformal time t — J a dt. Equation (140)) shows that the 
scalar degree of freedom has the effective mass 



1 _// 

M 2 = - — ^ = 

a 2 z s 



Q 2 S Qs ZHQs 

Ws'Ws^ 2Qs 

72F 2 H 4 



1 / 288# 3 - 12HR 1 
-F I r- + 



(2FH + f, RR R) 2 3 \2FH + f MR R f, RR/ 



+ (41) 



where we have eliminated the term H by using the background equations. 
In Fourier space the perturbation u satisfies the equation of motion 



/' + (fc 2 + Mia 2 ) 



0. 



(42) 



For k 2 /a 2 » M 2 , the propagation speed c s of the field u is equivalent to the speed of light c. Hence, in f(R) gravity, 
the gradient instability associated with negative c 2 is absent. For small k satisfying k 2 /a 2 <C M 2 , we require that 
M 2 > to avoid the tachyonic instability of perturbations. The viable dark energy models based on f(R) theories 
need to satisfy the condition RJ^rr -C F (i.e. m = Rf, RR /f,R <C 1) at early cosmological epochs in order to have 
successful cosmological evolution from radiation domination till matter domination. At these epochs the mass squared 
is approximately given by 



Mi 



■if. 



RR 



Under the no-ghost condition (l39l) the tachyonic instability is absent for 

f,RR > . 



(43) 



(44) 



The viable f(R) dark energy models have been constructed to satisfy the conditions (|39l) and (l4"4")l in the regime 
R > Ri, where R\ is the Ricci scalar at the late-time de Sitter point. Moreover we require that the models are 
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consistent with the conditions (1251) and (|27j) . The model (|28[) can be consistent with all these conditions, but the 
local gravity constraints demand that the variable m is very much smaller than 1 in the regions of high density (i.e. 
R Rq, where Rq is the cosmological Ricci scalar today). In the model (f28| one has m ~ n(— r — 1) around r ~ — 1. 
For the consistency with local gravity constraints we require that n < 10~ 10 [78j . but in this case the deviation from 
the ACDM model around the present epoch (R as Rq) is very small. 

If the variable m behaves as m = C(~r — l) p with p > 1 in the region R 3> i?o, then it is possible to satisfy local 
gravity constraints (i.e. m <C 0.01 ~ 0.1 for R 3> Rq) while at the same time showing deviations from the ACDM 
(m > 0.01 ~ 0.1 for R ~ i?o)- The models constructed in this vein are 



(A) f(R) = R- tiRc ,„,„ , 1 with ^ > °' " > and i? c > , (45) 



(fl/fl c ) 2 " 
' (R/R c ) 2n + 1 

(B) /(i?) = - (uEc [l - (l + R 2 /R 2 C ) with /i > 0, n > and R c > , (46) 



which were proposed by Hu and Sawicki [68( and Starobinsky [69(, respectively. R c is roughly of the order of the 
present cosmological Ricci scalar Rq for \x and n of the order of unity. The models (A) and (B) asymptotically behave 
as 

f(R) ~ E-//i? c [l - (R 2 /R 2 c )- n ] for i?>i? c , (47) 

which gives m(r) = C(— r — l) 2n+ . 

Another viable model that leads to the even rapid decrease of m toward the past is [Til] 

(C) f(R) =R- /ii? c tanh (R/R c ) with (i > 0, i? c > . (48) 

Other similar models were proposed by Appleby and Battye [70| and Linder [72j]. 
In what follows we shall discuss local gravity constraints on the above models. 



C. Local gravity constraints on f(R) gravity models 



Let us proceed to discuss local gravity constraints on f(R) gravity models. In the region of high density like Earth 
or Sun, the Ricci scalar R is much larger than the background cosmological value Rq . In this case the linear expansion 
of R = Rq + SR cannot be justified. In such a non-linear regime the chameleon mechanism (8~il.[8 ^| plavs an important 
role for the f(R) models to satisfy local gravity constraints [6 ll |64[ l68l l77H80l j (see also Refs. [184 -190]). 

To discuss the chameleon mechanism in f(R) gravity, it is convenient to transform the action ([1]) to the so-called 
Einstein frame action via the conformal transformation 11911 : 



_ q2 



9/ 



(49) 



The action in the Einstein frame includes a linear term in R, where the tilde represents quantities in th e Ei nstein 
frame. Introducing a new scalar field (j> = v /3/2«: 2 lnF, we obtain the action in the Einstein frame, as 0, Il9l| 



Se 



d 4 x 



J^R-lg^d^d^-V^) 



J tfxCuiF-^g^M), 



(50) 



where 



v(4>) 



RF-f 
2k 2 F 2 



(51) 



In the Einstein frame the scalar field <f> directly couples with non-relativistic matter. The strength of this coupling 
depends on the conformal factor f2 = y/F. We define the coupling Q as 



SI 2F V6' 



(52) 



which is of the order of unity in f(R) gravity. If the field potential V(<p) is absent, the field propagates freely with a 
large coupling Q. Since a potential (|51[) with a gravitational origin is present in f(R) gravity, it is possible for f(R) 
dark energy models to satisfy local gravity constraints through the chameleon mechanism [6ll HH, Wa , l80j . 
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In a spherically symmetric space-time under a weak gravitational background (i.e. neglecting the backreaction of 
gravitational potentials) , variation of the action (|50")l with respect to the scalar field (f> leads to 

d^ + 2d0 = dKff 
df 2 f df dtp ' 

where f is a distance from the center of symmetry, and V c s (</>) is an effective potential defined by 

Kff(0) = n0) + e Q V- (54) 

Here p* is a conserved quantity in the Einstein frame, which is related to the density p in the Jordan frame via the 
relation p* = e 3< ^^p. By the end of this section we use the unit n 2 = 1. 

We assume that a spherically symmetric body has a constant density p* — pa inside the body (f < r c ) and that 
the density outside the body (r > r c ) is p* = ps- The mass M c of the body and the gravitational potential $ c at 
the radius r c are given by M c = (4:7r/3)f^pA and <& c = M c /8nf c , respectively. The effective potential V e g((j)) has two 
minima at the field values 4>a and 4>b satisfying V^ s {4>a) — and V^. s {4>b) = 0, respectively (here a prime represents a 
derivative with respect to </>). The former corresponds to the region with a high density that gives rise to a large mass 
squared m A = V^ s (^>a): whereas the latter to the lower density region with a smaller mass squared m B = V^(</>b). 
When the "dynamics" of the field 4> with the field equation (|53|) is studied, we need to consider the effective potential 
(—V e &) so that it has two maxima at <f> = 4>a and <fi = 4>b- 

We impose the two boundary conditions (d<j>/df)(f = 0) = and cf>(r — > oo) — <\>b- The field (j> is at rest at f = 
and begins to roll down the potential when the matter-coupling term Qpac^ becomes important at a radius f\ in 
Eq. (j5"3")l . As long as fi is close to r c such that Af c = f c — fi <C f c , the body has a thin-shell inside the body. The field 
acquires a sufficient kinetic energy in the thin-shell regime (fi < r < f c ) and hence the field climbs up the potential 
hill outside the body (f > f c ). 

The field profile can be obtained by matching the solutions of Eq. ([53]) a t the radius f = r% and f — r c . Neglecting 
the mass term tub, the thin-shell field profile outside the body is given by jl92| 



where 



CM 

(f) - 4>b - 2Q eff ^ , (55) 



QeB ^ 3Qe t h , eth = — ■ (56) 

Here e t h is called a thin-shell para meter . Under the conditions Ar c /r c 1 and l/(mAr c ) -C 1, the thin-shell 
parameter is approximately given by |l92| 

£th ^ + — . (57) 
r c m A r c 

Provided that eth 1, the amplitude of the effective coupling Q c g can be much smaller than 1. It is then possible 
for the f(R) models (|Q| = l/\/6) to be consistent with local gravity experiments. Originally the thin-shell solution 
was derived by assuming that the field is frozen in the region < r < f i jsH, HH • In this case the thin-shell parameter 
is given by eth — Af c /f c , which is different from Eq. (|57[) . However, this difference is not important because the 
condition Ar c /r c l/(mAr c ) is satisfied for most of viable models [l92| . 

Consider the bound on the thin-shell parameter from the possible violation of equivalence principle (EP). The 
tightest bound comes from the solar system tests of weak EP using the free-fall acceleration of Moon (om ooh) and 
Earth (a®) toward Sun [82l |. The experimental bound on the difference of two accelerations is given by 193] 

|a M oon-ae| < 1Q -i3 (58) 



(«M + o®)/2 



Provided that Earth, Sun, and Moon have thin-shells, the field profiles outside the bodies are given by Eq. (f5"5"]) with 
the replacement of corresponding quantities. The acceleration induced by a fifth force with the field profile <p(r) and 
the effective coupling Q e g is a fifth = |Q c ffV</>(f)|. Using the thin-shell parameter e t h,© for Earth, the accelerations 
and ctMoon toward Sun (mass M ) are 



a© — ~o 



l + 18Q 2 e t 2 h , ffi |® 



GMr. 



"Moon — — 33 — 



1 + 18Q2e ^dir 



0*Moon 



(59) 
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7.0 x 10- 10 , and $ M oon 



where $ Q ~ 2.1 x 10 _b , $ e 
and Moon, respectively. Then the condition 



3.1 x 10 11 are the gravitational potentials of Sun, Earth 
translates to 



< 8.8 x HT 7 /|Q| 



(60) 



Since the condition \4>b\ \<t>A\ is satisfied for viable f(R) models (as we will see below), we have eth,® — 4>b /(6Q$®) 
from Eq. (|56|l. Hence the condition (jfJO)) corresponds to 



\4>i 



<3.7x 10~ 15 . 



(61) 



Let us consider local gravity constraints on the f(R) models given in Eqs. (|45l and (|46|) . In the region of high 
density where local gravity experiments are carried out, it is sufficient to use the asymptotic form given in Eq. (|47l) . 
In order for these models to be responsible for the present cosmic acceleration, R c is roughly the same order as the 
cosmological Ricci scalar Rq today for p and n of the order of unity. For the functional form f|47[) we have the following 
relations 



F 

v eS (4>) 



= 20/V6 



l-2np(R/R c )^ 2n+1 \ 



\pR c e~^^ 



1 - (2n + l) 



2n/(2n+l)' 



\\ffSnp J 

Inside and outside the body the effective potential (|63|) has minima at 

4>A * -VEnp(R c /p A ) 2n+1 , B ~ -V6np(R c /p B ) 2n+1 . 

If pa > pBt then one has \4>b\ > \<t>A\- 
The bound (IfJTT) translates into 



(62) 
(63) 

(64) 



rip 



n \ 2n+l 

^i) <1.5xHT 15 



(65) 



Here Xi is defined by x x = Ri/R c , where Ri is the Ricci scalar at the late-time de Sitter fixed point Pi given in 
Eq. (|18l) . Let us consider the model described by the Lagrangian density (l4"T|) for R> R\. If we use the models (|4"5")) 
and (|46p , then there are some modifications for the estimation of R\ . However this change is not significant when we 
place constraints on model parameters. 



The de Sitter solution for the model (j47|) satisfies p 
Eq. IBS]), it follows that 



2n+l 



/[2(xf 



2(x? 



Ri 

PB 



2n+l 



< 1.5 X 10 



-15 



1)]. Substituting this relation into 



(66) 



For the stability of the de Sitter point we require that m(Ri) < 1, which translates into the condition x\ n > In +3n+l 



Hence the term n/[2{x\ 



2/i 



1)] in Eq. (|66|) is smaller than 0.25 for n > 0. 



Let us use the simple approximation that Ri and ps are of the orders of the present cosmological density 10 -29 
g/cm 3 and the baryonic/dark matter density 10 -24 g/cm 3 in our galaxy, respectively. From Eq. (|66p we obtain the 
constraint I78l 



n > 0.9. 



(67) 



(1281) . it was shown in Ref. 
ACDM model is very small 



Thus n is not required to be much larger than unity. Under the condition (|67p . as R decreases to the order of R c , 
one can cosmologically see an appreciable deviation from the ACDM model. The deviation from the ACDM model 
appears when R decreases to the order of R c . The model (|48p also shows similar behavior. If we consider the model 

781 ] that the bound (|61l) gives the constraint n<3x 10~ 10 . Hence the deviation from the 
The models (|45|) and ((46]) are carefully constructed to satisfy local gravity constraints, 
while at the same time the deviation from the ACDM model appears even for n = 0(1). Note that the model (|4"5)l 
can easily satisfy local gravity constraints because of the rapid approach to the ACDM in the regim e R^S> R c . 

In the strong gravitational background (such as neutron stars), Kobayashi and Maeda |l94[ Il95| pointed out that 
for the model (|46[) it is difficult to obtain thin-shell solutions inside a spherically symmetric body with con stant 
density. For chameleon models with general couplings Q, a thin-shell field profile was analytically derived in Ref. [196| 
by employing a linear expansion in terms of the gravitational potential $ c at the surface of a compact object with 
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constant density. Using the boundary condition set by analytic solutions, Ref. 196] also numerically confirme d the 
existence of thin-shell solutions for $ c < 0.3 in the case of inverse power-law potentials V(4>) — Al 4+n (f)~ n . Ref. [l97l | 
also showed that static relativistic stars with constant density exists for the model (J46j) . The effect of the relativistic 
pressure is important around the center of the body, so that the field tends to roll down the potential quickly unless 
the boundary condition is carefully ch osen . Rea listic stars have densities pA (r) that globally decrease as a function of 
r. The numerical simulation of Refs. [l98l Il99j showed that thin-shell solutions are present for the f(R ) mo del (j46|) 
by considering a polytropic equation of state even in the strong gravitational background (see also Ref. [200J). 



III. SCALAR-TENSOR GRAVITY 



There is another class of modified gravity called scalar-tensor theories in which the Ricci scalar R is coupled to a 
scalar field ip. One of the simplest examples is the so-called Brans-Dicke theory with the action 

+ J d A xL M {g^^M), (68) 

where wbd is a constant (called the Brans-Dicke parameter), U(ip) is a field potential, and Cm is a matter Lagrangian 
that depends on the metric g^ v and matter fields \l/ m . The original Brans-Dicke theory [73j does not have the field 
potential. As we will see below, metric f(R) gravity discussed in Sec. [IT] is equivalent to the Brans-Dicke theory with 
w B d = 0. 



S = / d i x^/=g~ 



A. Scalar-tensor theories and the matter coupling in the Einstein frame 



The general action for scalar-tensor theories can be written as 



S 



d xy/^g 



l -f^R)- l -C{ V ){V V f 



J & A xC M (g^^M) , 



(69) 



the conformal factor fl = v F, the action (|69| can be transformed to that in the Einstein frame: 

1^ 1,.=. 



Se = 



-g 



iR-i{VcpY~v{4>) 



where / depends on the scalar field ip and the Ricci scalar R, ( is a function of <p. We choose the unit k 2 = SttG = 1. 
The action (|6"9"|) covers a wide variety of theories such as f(R) gravity (f(<p,R) — f(R), C — 0), Brans-Dicke theory 
(/ = ipR and ( = wbd/V), and dilaton gravity (/ = e~ v R and £ = — e~ v ). 
Let us consider theories of the type 

f(<p,R) = F(<p)R-2U(<p). (70) 
In order to avoid the appearance of ghosts we require that F(tp) > 0. Under the conformal transformation (|4Q|) with 



(71) 
(72) 

(73) 
(74) 



where 

V = U/F 2 . 

We have introduced a new scalar field <b in order to make the kinetic term canonical: 



dtp - 



'3 (F 



F 



We define the coupling between dark energy and non-relativistic matter in the Einstein frame: 



Q 



2F 



F, 



3 fF L 



F 



-1/2 



Recall that in metric f(R) gravity we have that Q = — l/v6- If Q is a constant, the following relations hold from 
Eqs. {Z3) and ((Til): 



F = e 



-2Q<f> 



( = (1-6Q')F[ ^ 



(75) 
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Then the action (1691) in the Jordan frame can be written as [76 



S = 



d 4 ^ 



-9 



F{4>)R 6Q 2 )F(</))(V0) 2 - Utf) 



J d 4 a; Cm{9^u, ^ m) 



(76) 



In the limit that Q — > 0, the action ([76)) reduces to the one for a minimally coupled scalar field <j) with the potential 
U (</)). The transformation of the Jordan frame action ([75)) via a conformal transformation g„ v = F{<f)g llv gives rise 
to the Einstein frame action (|71l) with a constant coupling Q. The action (fTTj) is equivalent to the action ([50)) with 

2Q</> 



One can compare (|76p with the action ([68| in Brans-Dicke theory. Scttin. 
actions are equivalent if the parameter wbd is related to Q via the relation [7j 

3 + 2w B d 



F 



2Q 2 



one finds that two 



(77) 



Using this relation, we find that the General Relativistic limit (wbd — > oo) corresponds to the vanishing coupling 
{Q — >• 0). Since Q = — in metric f(R) gravity, this corresponds to the Brans-Di cke param eter ojbd = [5lL Wty. 
The experimental bound on cjbd for a massless scalar field is given by ojbd > 40000 [193, 201], which translates into 
the condition 



|Q| < 2.5 x 10" 



(for a massless field). 



(78) 



In such cases it is difficult to find a large difference relative to the uncoupled quintessence model. In the presence of 
the field potential, however, it is possible for large coupling models (|Q| ~ 1) to satisfy local gravity constraints via 
the chameleon mechanism [76j. 

The above Brans-Dicke theory is one of the examples in scalar-tensor theories. In general the coupling Q is field- 
dependent apart from Brans-Dicke theory. If we consider a nonminimally coupled scalar field with F(ip) = 1 — ^ip' 2 
and (((p) = 1, then it follows that Q(<f>) = £<ft /[l — 1 — 6£)1 1// 2 . The cosmological dynamics in such a theory have 
been studied by a number of authors [H-lifl 120^42051 l207M2lTI |. If the field is nearly massless during most of the 



cosmological epochs, the coupling Q needs to be suppressed to avoid the propagation of the fifth force. 

In the following we shall study the cosmological dynamics and local gravity constraints on the constant coupling 
models based on the action ([76i) with F{4>) = e~ 2Q *. 



B. Cosmological dynamics in Brans-Dicke theory 



We study the cosmological dynamics for the Jordan frame action (|76j) in the presence of a non-relativistic fluid with 
energy density p m and a radiation fluid with energy density p r . We regard the Jordan frame as a physical frame due 
to the usual conservation of non-relativistic matter (p m cx a -3 ). In the flat FLRW background variation of the action 
([TB")) with respect to g^ v and 4> gives the following equations of motion 



3FH 2 = (1 - 6Q 2 )F4> 2 /2 + U- 3HF + p m + p r , 
2FH = -(1 - 6Q 2 )F</> 2 - F + HF - p m - (4/3)p r 



(1-6Q 2 )F[ 
Let us introduce the following variables 



3H<j> + F/(2F)4>] + + QFR = . 



Xi 



and 



ZFH 2 ' 

These satisfy the relation VL m + f2 ra d 
equations for xi, X2 and x$: 



^rad 
--1)1 



— — 

H\3F 



X3 



1 

H 




Pr 

3F ' 



(79) 
(80) 
(81) 

(82) 



n BE = {l-6Q 2 )x 2 +x 2 +2V6Qx 1 . (83) 
1 from Eq. ([79]). Using Eqs. ([79 ]) - ([81 ) . we obtain the differential 



dxi y/6 2 fn \ V6Q 



(5 - 6Q 2 )xl + 2V6Q Xl - 3a 



1 



Xi 



H 

FT 2 



dx 2 
dN 
dx 3 
dN 



V6 



(2Q - \)xix 2 - x 2 



V&Qxix 3 - 2x 3 



%3 



H 

m 

H 



(84) 
(85) 
(86) 
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where TV — In a, A — —U^/U, and 



H_ 1 -6Q 2 
TP- ~ 2 



3 + 3x 2 - 3a; 2 , + x 2 - 6Q 2 x 2 + 2V6Qxi) + ZQ(\x\ - AQ) . 



(87) 



The effective equation of state of the system is given by by w c s = — 1 — 2H / (3H 2 ). 

If A is a constant, one can derive the fixed points of the system (|84 |) -(|86 p in the absence of radiation (2:3 = 0) [76| : 



(a) 



VGQ A 3-2Q 2 _ 4Q 2 

3(2Q 2 -1)' )' 3(1 -2Q 2 ) 2 ' Wcff 3(1 -2Q 2 ) 



(88) 



(b) 



[xi,x 2 ) = [ —;=—— — ,0 ] , S2 m = 0, w eff - 



VV6Q±1 : 



3(1 ± V6Q) 



(89) 



(c) 



(an, a?2) 



76(4Q - A) 



6(4Q 2 -QA- 1)' [ 6(4Q 2 -QA-1 



6 - A 2 + 8QA - 16Q 2 



(d) 



1/2N 



n m = o , tDcff = - 



20Q 2 - 9QA - 3 + A 2 
3(4Q2-QA-1) 



(90) 



, , f V6 /3 + 2QA-6Q 2 \ 1 3 - 12Q 2 + 7QA 2Q 
(^^)=Ua'V 2^ ' fim = 1 V ' ^ = " — ■ 



(91) 



(e) 



(xi,x 2 ) = (0,1) , f2 m = 0, w c ff = -1 



(92) 



The point (e) corresponds to the de Sitter point, which exists only for A = AQ [this can be confirmed by setting — 
in Eqs. CSJ-dST)]. 

We first study the case of non-zero values of Q with constant A, i.e. for the exponential potential U(cf>) — Uoe~ x ^. 
We do not consider the special case of A = AQ. The matter-dominated era can be realized either by the point (a) or 
by the point (d) . If the point (a) is responsible for the matter era, the condition Q 2 <C 1 is required. We then have 
fi m ~ 1 + 10Q 2 /3 > 1 and w c s — 4Q 2 /3. When Q 2 <C 1 the scalar-field dominated point (c) yields an accelerated 
expansion of the Universe provided that — v2 + AQ < A < v2 + AQ. Under these conditions the point (a) is followed 
by the late-time cosmic acceleration. The scaling solution (d) can give rise to the equation of state, w c s — for 
\Q\ <C |A|. In this case, however, the condition w c g < —1/3 for the point (c) gives A 2 < 2. Then the energy fraction of 
the pressureless matter for the point (d) does not satisfy the condition O m ~ 1. From the above discussion the viable 
cosmological trajectory for constant A corresponds to the sequence from the point (a) to the scalar-field dominated 
point (c) under the conditions Q 2 <C 1 and — spl + AQ < A < \/2 + AQ. 

We shall proceed to the case where A varies with time. The fixed points derived above for constant A can be 
regarded as the "instantaneous" fixed points, provided that the time scale of the variation of A is smaller than that of 
the cosmic expansion. The matter era can be realized by the point (d) with \Q\ <C |A|. The solutions finally approach 
either the de Sitter point (e) with A = AQ or the accelerated point (c). 

In the following we focus on the case in which the matter era with the point (d) is followed by the accelerated 
epoch with the de Sitter solution (e). To study the stability of the point (e) we define a variable x± = F, satisfying 
the following equation 



dx4 
dN 



= -2V6Qxi 



X4 . 



(93) 
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Considering the 3x3 matrix for perturbations Sxi, 6x2 and Sx^ around the point (e), we obtain the eigenvalues 

-3, - 



dA 



l±^jl-- Fl Q— (Fr) 



(94) 



where F± = F(4>i) is the value of F at the de Sitter point with the field value 0i. Since F\ > 0, we find that the de 
Sitter point is stable under the condition 



Q^(Fi)>0, i.e. 



dA 



l)<0. 



(95) 



Let us consider the f(R) model (j47|) in which the models (|45|) and (j46|) are recovered in the regime R^> R c . Since 
e 24>/Ve = 1 _ 2n M (i?/i? c )-(2«+ 1 ), the potential U = (FR - f)/2 is given by 



£/(0) 



2n + 1 



(2npJ 



±i u 

n/(2ri+l) ^ 



2<j>/s/E 



2n/(2n+l) 



In this case the slope of the potential, A = —U^/U, is 



A 



4ne 20/V6 



%/6(2n^) 2 "/(2n 



+1) 



2n + 1 



(2n/i) ; 



n/(2n+l) V 



-1 -2n/(2n+l) 



(l - e 2< ^) 



-l/(2n+l) 



(96) 



(97) 



In the deep matter-dominated epoch during which the condition R/R c 3> 1 is satisfied, the field is very close to 
zero. For n and p of the order of unity, we have |A| ^> 1 at this stage. Hence the matter era can be realized by the 
instantaneous fixed point (d). As R/R c gets smaller, |A| decreases to the order of unity. If the solutions reach the 
point A = AQ = — 4/V6 and satisfy the stability condition dA/dF < 0, then the final attractor corresponds to the de 
Sitter fixed point (e). 

For the theories with general couplings Q, it is possible to construct a scalar-field potential that is the generalization 
of (|5Tj) . One example is (7(| 



U{4>) = U [1- C(l - 



-2Q<j>\p 



(U >0, C > 0, < p < 1) . 



The f(R) model ([47]) corresponds to Q = — l/v6 and p = 2n/(2n + 1). The slope of the potential is given by 

2CpQe- 2 W(l-e- 2 W)f- 1 



A 



1 - C(l - e~ 2 W)P 



(98) 



(99) 



We have U{4>) -> U for -> and C/(0) [/ (1 - C) in the limits -> 00 (for Q > 0) and -00 (for Q < 0). 

The field is nearly frozen around the value = during the deep radiation and matter epochs. In these epochs we 
have R ~ p m /F from Eqs. (fT9"|) - ([8"Tj) by noting that C/q is negligibly small compared to p m or p r . Using Eq. (I5T1) . it 
follows that f7,£ + Qpm ~ 0. Hence, in the high-curvature region, the field evolves along the instantaneous minima 
given by 



J_ ( 2U oP C \ 

2Q V Pm ) 



V(i-p) 



(100) 



The field value |0 m | increases for decreasing p m . As long as the condition p m 3> 21] ^pC is satisfied, we have |0„,.| <C 1 
from Eq. (fTOO]) . 

For field values around = one has |A| S> 1 from Eq. (|99|). Hence the instantaneous fixed point (d) can be 
responsible for the matter-dominated epoch provided that \Q\ <C |A|. The variable F = e~ 2 ^ decreases in time 
irrespective of the sign of the coupling Q and hence < F < 1. The de Sitter solution corresponds to A = AQ, that is 



C 



(1 - FxJp-i [2 + (p - 2)F X ] 



(101) 



This solution is present as long as the solution of this equation exists in the region < F\ < 1 . 
From Eq. (|99|) the derivative of A with respect to is 

dA _ 4C P Q 2 F(1 - F)p- 2 [1 - P F- C(l - F) p ] 
d0 ~ [1 -C(l -F)p] 2 ' 



(102) 
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Figure 2: The evolution of SIde, f2 ra( j and w c g in Brans-Dicke theory with the potential (|98[) . The model parameters are 
Q = 0.01, p = 0.2 and C* = 0.7 with the initial conditions x\ = 0, x 2 = 2.27 x 10~ 7 , x 3 = 0.7, and a: 4 - 1 = -5.0 x 10~ 13 . From 
Ref. [H. 



The de Sitter point is stable under the condition 1 — pF\ > (7(1 — -Fi) p - Using Eq. (I101[) this translates into 

F x > 1/(2 -p). (103) 

When < C < 1 one can show that d\/d<f) < is always satisfied. Hence the solutions approach the de Sitter 
attractor after the end of the matter era. When C > 1, the dc Sitter point is stable under the condition (|103[) . If this 
condition is violated, the solutions choose another stable fixed point [such as the point (c)] as an attractor. 

The above discussion shows that, if < C < 1, the matter point (d) can be followed by the stable de Sitter solution 
(e). In Fig. [2] we plot the evolution of JIde, ^m, ^rad, and w e g for Q = 0.01, p = 0.2 and C = 0.7. This shows that 
the viable cosmological trajectory can be realized for the potential (f9"5)h In order to confront with SN la observations, 
it is possible to rewrite Eqs. (|T5|) and (15U1) in the forms of Eqs. (|29p and (|3H)l by defining the dark energy density 
Pde and the pressure Pde hi the similar way. It was shown in Ref. [76j that the phantom equation of state as well 
as the cosmological constant boundary crossing can be realized for the field potentials U(c/>) satisfying local gravity 
constraints. 



C. Local gravity constraints on Brans-Dicke theory 

We study local gravity constraints on Brans-Dicke theory described by the action ([7d| . In the absence of the 
potential U(cf>) we already mentioned that the Brans-Dicke parameter wbd is constrained to be wbd > 4.0 x 10 4 from 
solar-system experiments. This gives the upper bound (|78p on the coupling Q between the field (f> and non-relativistic 
matter in the Einstein frame. This bound also applies to the case of a nearly massless field with the potential U(cj>) in 
which the Yukawa correction e~ Mr is close to unity (where M is the scalar field mass and r is an interaction length). 

In the presence of the field-potential it is possible for large coupling models (|Q| ~ 1) to satisfy local gravity 
constraints provided that the mass M of the field <j> is sufficiently large in the region of high density. In fact, 
the potential ([55]) is designed to have a large mass in the high-density region, so that it can be compatible with 
experimental tests of gravity through the chameleon mechanism. In the following we study the model (I98[) and derive 
the conditions under which local gravity constraints can be satisfied. If we make a conformal transformation for the 
action (|98|) . the action in the Einstein frame is given by ([76| with F(<f>) — e~ 2 ^ . We can use the results obtained in 
Sec. Ill C[ because thin-shell solutions have been derived for the general coupling Q. 
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As in the case of f(R) gravity, we consider a configuration in which a spherically symmetric body has a constant 
density pA inside the body and that the density outside the body is given by p — ps (^C Pa)- Under the condition 
\Q4>\ <C 1, we have ~ —2UoQpC{2Q(j)) p ^ 1 for the potential V — U/F 2 in the Einstein frame. Then the field values 
at the potential minima inside and outside the body are 

, if^y ,,, - rt , , i»y /<, -'\ (104) 



2Q V PA J ' r 2Q V PB 

In order to realize the accelerated expansion today, the energy scale Uq is required to be the same order as the square 
of the present Hubble parameter Ho, i.e. Uq ~ Hq ~ po, where po ~ 1CP 29 g/cm 3 is the cosmological density today. 
The baryonic/dark matter density in our galaxy corresponds to ps — 1CP 24 g/cm 3 . Hence the conditions IQ^aI <C 1 
and |<30b| <C 1 are in fact satisfied unless C ^> 1. The field mass squared m 2 A = at = <f>A is approximately 
given by 



<-w^m^ Q 'ii) u °- (105) 

which means that tua can be much larger than Ho because of the condition pa Uo- This large mass allows the 
chameleon mechanism to work, because the condition 1/(771^4^) -C 1 is satisfied. 

The bound (|61|) coming from the violation of equivalence principle in the solar system translates into 

(2UopC/p B ) 1/{1 - p) < 7.4 x 1CT 15 \Q\ . (106) 

Let us consider the case in which the solutions finally approach the dc Sitter solution (c). At the point (e), one has 
3FxHf = U [l - C(l - Fi) p ] with C given in Eq. ([TDT]) . Hence we get the following relation 

U = 3H*[2 + (p-2)F 1 ]/p. (107) 

Plugging this into Eq. (|106|) . it follows that 

(i?i/PB) 1/(1 " p) (1 - F 1 ) < 7.4 x 10 15 |Q| , (108) 

where R\ = \2H\ is the Ricci scalar at the dc Sitter point. Since the term (1 — F\) is smaller than 1/2 from the 
condition ((T03|) . we obtain the inequality (Ri/ p B ) 1/{1 ~ p) < 1-5 x 10~ 14 |<2|. Using the values i?i = 10~ 29 g/cm 3 and 
Pb = 10 -24 g/cm 3 , we obtain the following bound 

^'- 13.8-Lior <109) 

When |Q| = 10 _1 and \Q\ = 1 we have p > 0.66 and p > 0.64, respectively. Thus the model can be compatible with 
local gravity experiments even for \Q\ = 0(1). 

In Ref. [79| it was shown that in order to satisfy both local gravity and cosmological constraints the chameleon 
potentials in the Einstein frame need to be of the form V((f>) = M 4 [l + f((f>)}, where the function f((p) is smaller 
than 1 today and M is a m ass that corresponds to the dark energy scale (M ~ 10~ 12 GeV). The potential V((j>) = 
M 4 exp[n(M / cf>) n ] [82|, [212] is one of those via ble candidates, but the allowed model parameter space is severely 
constrained by the 2006 Eot-Wash experiment [213J. Unless the parameter p is unnaturally small (p < 10 -5 ), this 
potential is incompatible with local gravity constraints for {n, Q} — 0(1). 

On the other hand, the chameleon potential V(4>) — Vo[l — p(l — e~^) n ] (0 < n < 1) can satisfy both local gravity 
and cosmological constraints 4 . In Ref. [79[ this potential is consistent with the constraint comi ng from 2006 Eot- 
Wash experiments as well as the WMAP bound on the variation of the field-dependent mass [214j for natural model 
parameters. 



4 In the Einstein frame the potential ((98j takes the form V(cf>) = U/F 2 = U e A ^^ [l - C*(l - e'^^Y] , so in the region \Q<f>\ < 1 this 
potential is similar to V{(j>) = Vo[l - fi(l - e~^)"]. 
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IV. DGP MODEL 



In this section we review braneworld models of dark energy motivated by string theory. In braneworlds standard 
model par t icles are confined on a 3-dimensional (3D) brane embedded in 5-dimensional bulk with large extra dimen- 
sions [215l |216| |. Dvali, Gabadadze, and Porrati (DGP) [4l[ proposed a braneworld model in which the 3-brane is 
embedded in a Minkowski bulk with infinitely large extra dimensions. On e ca n recover Newton's law by adding a 4D 
Einstein-Hilbcrt action sourced by the brane curvature to the 5D action |217l |. The presence of such a 4D term may 
be induced by quantum corrections coming from the bulk gravity and its coupling with matter on the brane. In the 
DGP model the standard 4D gravity is recovered at small distances, whereas the effect from the 5D gravity manifests 
itself for lar ge d i stan ces. Interesting ly on e can realize the self cosmic acceleration without introducing a dark energy 
component [2ll|21jj (see also Ref. (220j |). 



A. Self-accelerating solution 



The action of the DGP model is given by 

1 



., . d 6 X y /-gR+— s - I d 4 X^R- I d b Xy/-gC M , (110) 

ZK (5) J ZK (4) 



where c/ab is the metric in the 5D bulk and g^ u — d ll X A d ll X B gAB is the induced metric on the brane with X A (x c ) 
being the coordinates of an event on the brane labelled by x c . The 5D and 4D gravitational constants, k? 5 n and t? 4 \, 
are related with the 5D and 4D Planck masses, Mtg\ and Mu-i, via 



2 -1/Mf 5)) 4) = 1/Mf 4) . (Ill) 



(5) 



The first and second terms in Eq. (jllOp correspond to Einstein-Hilbert actions in the 5D bulk and on the brane, 
respectively. There is no contribution to the Lagrangian Cm from the bulk because we are considering a Minkowski 
bulk. Then the matter action consists of a brane-localized matter whose action is given by J d A Xyf—g (<r + £^ ano ), 
where a is the 3-brane tension and C M anc is the Lagrangian density on the brane. Since the tension is unrelated to 
the Ricci scalar R, it can be adjusted to be zero (as we do in the following). 

In order to study the cosmological dynamics on the brane (located at y — 0), we take a metric of the form: 

ds 2 = -n 2 {T,y)dr 2 + a 2 (r,y) 'y ij dx i dx j + dy 2 , (112) 

where jij represents a maximally symmetric space-time with a constant curvature K. The 5D Einstein equations are 

Gab = Rab - ^R§ab = Kf 5 )TAB , (113) 

where Rab is the 5D Ricci tensor, Tab is the sum of the energy momentum tensor jj^ ane ) on the brane and the 
contribution XJ ab coming from the scalar curvature of the brane: 

Tab = T ( A h r c) + V ab ■ (114) 
Since we are considering a homogeneous and isotropic Universe on the brane, one can write 7 1 ^ branG '' m the form 

T# (branc) - 8(y) diag(- PM , P M , Pm,Pm,0) , (115) 

where pu and Pm are functions of r only. The non-vanishing components coming from the Ricci scalar R of the 
brane are 

U 00 = -4-(^+K^) S(y), (116) 
kt.n \a z a z I 



"(4) 
K (4 ) 



— I ^ + 2 2-| A 

a z an a 



m s(v) , (H7) 
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whe re a dot repre sents a derivative with respect to r. The non- vanishing components of the 5D Einstein tensor Gab 
are HTl EHl SH 



Goo 
dj - 

G55 



a n 
2— + — 

a n 



. an a 

an a 



a' 2 N 


n 2 ~ 

+ ^ 




a 2 y 




a' 2 


an J 


a 2 


— H 


+ — 

n z 



ii a' 2 
-2 - 

1- 



an 
an 



a 



3^ 

2 



an 
2 — 

an 



A" 



(118) 
(119) 
(120) 
(121) 



where a prime represents a derivative with respect to y. 

Assuming no flow of matter along the 5-th dimensions, we have T05 = and hence G05 = 0. Then Eqs. (|118p and 
(|121[) can be written as 



where 



3n 2 

Goa = ~2M 1 '' 



(a'a) 2 



G55 — — 



2a 3 d 



(da) 2 



Ka 1 



(122) 



(123) 



Since we are considering the Minkowski bulk, we have Goo 
1 = 0. Integrations of these equations lead to 



( a ' a f _ l^L _ Ka 2 
n z 



and G55 = locally in the bulk. This gives I' = and 
G = 0, (124) 



where C is a constant independent of r and y. 

We shall find solutions of the Einstein equations (1113j) in the vicinity of y = 0. The metric needs to be continuous 
across the brane in order to have a well-defined geometry. However, its derivatives with respect to y can be discon- 
tinuous at y — 0. The Einstein tensor is made of the metric up to the second deriva tives with respect to y, so the 
Einstein equations with a distributional source are written in the form [218l l22ll . |222| 

(125) 



g" = TS(y), 

where 6(y) is a Dirac's delta function. Integrating this equation across the brane gives 

[g']=T, where [g'\ = g'(0+) - «/((T) . 



(126) 



The jump of the first derivative of the metric is equivalent to the energy-momentum tensor on the brane. 

Equations (|118[) and (|1 19[) include the derivatives a" and n" of the metric. Integrating the Einstein equations 



Goo — K (5)^oo an d Gj 



K (5)-^j across the brane, we obtain 



a b 

K] 

n b 



K 2 K 2 

K (5) n , K (5) 



— + K— 

a b a b 



; Ofc nb_ 
ab nb 



,a& 
a h 



K- 



(127) 
(128) 



where the subscript "6" represents the quantities on the brane. 

We assume the symmetry y -f-> — y, in which case [a'] = 2a'(0 + ) and [n'] = 2n'(0 + ). Substituting Eq. (|127|) into 
Eq. (|124j) . we obtain the modified Friedmann equation on the brane: 




^(5) 



K 

72 



Pa/ ■ 



(129) 
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where H = ^/(atnt) is the Hubble parameter a nd e = ±1 is th e sign of [a 1 ]. The constant C can be interpreted as 
the term coming from the 5D bulk Weyl tensor [218, 2191 1223| . Since the Weyl tensor vanishes for the Minkowski 
bulk, we set C = in the following discussion. We introduce a length scale 



r c = tt-^t - = ttttt - ■ (130) 

2 «(4) 2M ( 5) 



Then Eq. (jl 29[) can be written as 



.2 



r c V a a o 



H* + — = H 2 + — - -^ PM , (131) 



where we have omitted the subscript "6" for the quantities at y — 0. 

Plugging the junction conditions (|127l) and (|128l) into the (05) component of the Einstein equations, G05 = 0, the 
following matter conservation equation holds on the brane: 

^- + 3H(p M + P M ) = 0, (132) 
at 

where t is the cosmic time related to the time r via the relation At = nbdr. If the equation of state, wm = Pm/pm, 
is specified, the cosmological evolution is known by solving Eqs. (|131[) and (|132[) . 
For a flat geometry (K — 0), Eq. (|131|) reduces to 

H*-±H=^p M . (133) 

If the crossover scale r c is much larger than the Hubble radius if -1 , the first term in Eq. (|133|) dominates over the 
second one. In this case the standard Friedmann equation, H 2 = k? 4 )Pm/3, is recovered. On the other hand, in the 

regime r c < H , the presence of the second term in Eq. (| 133|) leads to a modification to the standard Friedmann 
equation. In the Universe dominated by non-relativistic matter (pm oc a~ 3 ), the Universe approaches a de Sitter 
solution for the branch e = +1: 

H^H dS = l/r c . (134) 
We can realize the cosmic acceleration today provided that r c is of the order of the present Hubble radius H^ 1 . 



B. Observational constraints on the DGP model and other aspects of the model 
Equation (I131[) can be written as 

K 



H 1 




1 



1 



3 PM + ^ + 2^. 



(135) 



For the matter on the brane, we consider non-relativistic matter with the energy density p m and the equation of state 
w m = 0. We then have p m = pm\l + z) 3 from Eq. (|132[) . Let us introduce the following density parameters 



n 



(0) 
K 



K 



n 2 H 2 



1 



ArlH 2 



K 2 O {0) 



Then Eq. (fT35]) reads 



H 2 (z) = H 2 



n { ° ] (1 + zf + ( y] fiW (1 + z)3 + + ^ 



(0) 



(136) 



(137) 



The normalization condition at z = is given by 

nw + nW + 2 V /i-nWv^ = 1 - 



(138) 
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Figure 3: Observational constraints on the DGP mo del f rom the SNLS data [2251 ] (solid thin), the BAO [f| (dotted), and the 
CMB shift parameter from the WMAP 3-year data [2261 ] (dot-dashed). The thick line represents the curve (|139f) for the flat 
model (n^ = 0). The fi gure labels Q, m and f2 rc correspond to Qm and , respectively. From Ref. [9lJ . 



For the fiat universe (K = 0) this relation corresponds to 

nCf? = (i-n(?) a /4. (139) 

The parametrization ([1371) of the Hubble parameter together with the normalization ([13811 can be used to place 
observational constraints on the DGP model at the background level [89T - I931 ] . In Ref. [89j| the authors found a 
significantly worse fit to Sup ernova la (SN la) data and the distance to the last-scattering surface [H relative to the 
ACDM model. In Refs. [90 and [HJ the authors showed that the flat DGP model is disfavored from the combined 
data analysis of SN la [224 . |225| and BAO @. In Fig. [3[w e sho w the joint observational constraints [9lJ from the 
data of SNLS [225j . BAO [|, and the CMB shift parameter [HI]. While the flat DGP model can be consistent with 
the SN la data, it is under strong observational pressure by adding the data of BAO and the CMB shift parameter. 
The open DGP model gives a slightly better fit relative to the flat model (9j |93|. The joint analysis using the data of 
SN la, BAO, CMB, gamma ray bursts, and the linear growth factor of matter perturbations show that the flat DGP 
model is incompatible with current observations (94{. 

In the DGP model a brane-bending mode <f> (i.e. longitudinal graviton) gives rise to a field self-interaction of the 
form □</>(<9 M </><9 M 0) through a mixing with the transverse graviton [83|,|84|. This can lead to the decoupling of the field 
cj> from gravitational dynamics in the local region by the so-called Vainshtein mechanism [85[ . The General Relativistic 
behavior can be recovered within a radius r* = (rgr^) 1 / 3 , where r g is the Schwarzschild radius of a sou rce. Since 
r* is larger than the solar-system scales, the DGP model can evade local gravity constraints [H, HH, |228 [. However 
the DGP model is plagued by a strong coupling problem for typical distances smaller than 1000 km [229j ]. Some 
regularization metho ds have b een proposed to avoid the strong coupling problem, such as smoothing out the delta 
profile on the brane [230l . l23l| or re-using the delta function profile but in a higher-dimensional brane [232ll234T| . 



As we will see in I VII C[ the analysis of 5D cosmological perturbations on the scales larger than r* shows that the 
DGP model contains a ghost mode in the scalar sector of the gravitational field [8614881 ]. There are several ways of 
the generalization of the DGP mo del t o avoid the appearance of ghosts. One way is to consider the 6D braneworld 



set-up as in the Cascading gravi ty 235j . Another is to generalize the field self-interaction term □0(<9 p </x9 A ,0) to more 



general forms in the 4D gravity [42[. In Sec. [V] we shall discuss the latter approach ("Galileon gravity") in detail. 
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V. GALILEON GRAVITY 



In the DGP model the field derivative self-coupling Dcf>(d ^ipd^' <f>) , arising from a brane-bending mode, allows the 
decoupling of the field from matter within a Vainshtein radius. In the local regions where solar-system experiments 
are carried out, the field is nearly frozen through the non-linear self-interaction. This is different from the chameleon 
mechanism in which the presence of the field potential with a matter coupling gives rise to a minimum with a large 
mass in the regions of high density. 

Under the Galilean shift d^cj) —> d^ip+b^, the field equation following from the Lagrangian D0(9^ 4>d '^(f>) is unchanged 
in the Minkowski space-time. The generalization of the nonlinear field Lagrangian to more general cases may be useful, 
e.g., to overcome the ghost problem associated with the DGP model. In fact Nicolis et al. [42| derived five Lagrangians 
that lead to the field equations invariant under the Galilean shift d^cf) — > d^4> + in the Minkowski space-time. The 
scalar field respecting the Galilean symmetry is dubbed "Galileon" . Each of the five terms only leads to second-order 
differential equations, keeping the th eory free from unstable spin-2 ghost degrees of freedom. 

If we extend the analysis in Ref. |42|| to that in the curved space-time, the Lagrangians should be promoted to 
covariant forms. Deffayet et al. [95l.l96j| derived covariant Lagrangians £; (i — 1, • • • ,5) that keep the field equations 
up to second-order, while recovering the five Lagrangians derived by Nicolis et al. in the Minkowski space-time. This 
can be achieved by introducing field-derivative couplings with the Ricci scalar R and the Einstein tensor G„ p in the 
expression of £4,5. Since the existence of those terms affects the effective gravitational coupling, the Galileon gravity 
based on the covariant Lagrangians Ci (i = 1, • • • ,5) can be classified as one of modified gravitational theories. 

The cosmological dynamics including the terms up to £4 and £5 have been studied by a number of authors [97l . l236j |. 
In particular Refs. (97|, [98| have shown that, for the covariant Galileon theory having de Sitter attractors, cosmological 
solutions with different initial conditions converge to a common trajectory- a tracker solution. Moreover there is a 
viable parameter space in which the conditions for the avoidance of ghosts and Laplacian instabilities of scalar and 
tensor perturbations are satisfied. 

The generalization of Galileon gravity, whi ch m o stly c orresponds to the modification of the term £3 = D^d^ipd^cf)), 
has been also extensively studied recently [237]- [258j |. One application is to introduce the non- li near field self - 
interaction of the form £ i ((/))D(p(d l j l (j)d tJ, (j)) in the action of (generalized) Brans-Dicke theories [l83l [238-240, |242| |. 
where £ is a function of <\>. For suitable choices of the function £(<^>), there exist de Sitter (dS) solutions responsible 
for dark energy even in the absence of the field potential. The cosmol ogy based on a further gener al term G(4>, X)Ocj) 
has been discussed in the context of either dark energy and inflation [246l [24 ^ [252 L [253 l l256H258j . 

In the following we review the cosmological dynamics in the Galileon dark energy model based on the covariant 
Lagrangians Ci (i — 1, • • • ,5) and study observational constraints on the model. We will also discuss the modified 
version of Galileon gravity in which the term £3 = D(j>(d fJl (/)d ll (f>) is generalized. 



A. Cosmology of a covariant Galileon field 



We start with the covariant Galileon gravity described by the action [95|, |96 



S 



d x^g 



R 



2 ^ 

i=l 



C l C l 



d 4 xC 



M ■ 



(140) 



where M p i = (8nG) x l 2 is the reduced Planck mass, and q's are constants. The five Lagrangians d (i = 1, • • • , 5) 
satisfying the Galilean symmetry in the limit of the Minkowski space-time are given by 

A = M 3 , £ 2 - (V0) 2 , £3 = (n0)(V0) 2 /M 3 , £4 = (V</>) 2 [2(D0) 2 - - #(V0) 2 /2] /M 6 , 

£ 5 = (V^) 2 [(D0) 3 - 3(D0) 0;^^ + 20 ; /0 ; /0 ; / - e^^G.pJ/M 9 , (141) 

where M is a constant having a dimension of mass, and G vp is the Einstein tensor. For the matter Lag rangian £m we 
take into account perfect fluids of non-relativistic matter (energy density p m , equation of state w m — 0) and radiation 
(energy density p r , equation of state w r — 1/3). 

Let us consider the FLRW metric with the cosmic curvature K: 

<l? "~ 2 (d0 2 +sin 2 #d0 2 ) . 



ds 2 = -dt z + a 2 (t) 



1 - Kr 2 

Variation of the action (|140p with respect to g^ v leads to the following equations of motion 

3M 2 Y H 2 = p DE + p m + Pr + PK , 

3M^H 2 + 2A/ 2 1 # 



-Pi 



DE 



Pr/3 + Pk/3 , 



(142) 

(143) 
(144) 
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where px = —SKM^/a 2 , and 

Pde = ~C!M 3 (l)/2-C20 2 /2 + 3c3H<j) 3 /M 3 -45c4H 2 <p 4 /(2M 6 ) + 2lc 5 H 3 <j) 5 /M 9 , (145) 
P DE = Cl M 3 (j)/2- c 2 2 /2- c^ 2 ip/M 3 + 3c 4 </> 3 [8ff0+ (3H 2 + 2i2>]/(2M 6 ) 

~3c 5 H<j> 4 [5H4> + 2{H 2 + H)<j)]/M 9 . (146) 

The matter fluids satisfy the continuity equations p rn + 3Hp m = and p r + 4Hp r = 0. We define the dark energy 
equation of state wde and the effective equation of state uu e s, as 

»DE= , WcS = -1 - ■ (147) 

PDE 3iP 

Using the continuity equation pde + 3-ff (pde + -Pde) = 0, it follows that wde = w e s — ^de/ (3-H^de)- 

Since we are interested in the case where the late-time cosmic acceleration is realized by the field kinetic energy, 
we set ci = in the following discussion 5 . Then the de Sitter solution (H — H d s = constant) can be present for 
4> = <AdS = constant. We normalize the mass M to be M 3 = M p iH% s , which gives M m KT 40 M p i for H ds m 10 _60 M p i. 
Defining x d $ = <fids/(HdsM p \), Eqs. (| 143|) and (|144[) lead to the following relations at the de Sitter solution: 

c 2 x% = 6 + 9a - 12(3 , c 3 :r ds = 2 + 9a - 9/3 , (148) 

where 

a = c 4 xj s , (3 = c 5 x 5 dS . (149) 

It is convenient to use the variables a and /3, because the coefficients of physical quantities and dynamical equations 
can be expressed by a and j3. The relations (|148[) do not change under the rescaling Xds — > "fx d s and Ci — > Ci/rf, 
where 7 is a real constant. Then the rescaled choice of Ci can provide the same physics. 

In order to study the cosmological dynamics, we introduce the following dimensionless variables: 

_ 0dS#dS _ 1 ( <j> \ , 1KrV . 

n = j— , r 2 = — H— • (150) 

At the de Sitter solution r% — 1 and r 2 = 1. We define the dark energy density parameter 

fi DE = = -(2 + 3a - 4/3)r?r 2 /2 + (2 + 9a - 9/3)r 2 r 2 - 15ar x r 2 /2 + 7/3r 2 . (151) 

oiW pl ii 

Then Eq. (|T43)) can be written as il DE + OT + Q r + Q K = 1, where Q ro = p m /(3M 2 y H 2 ) : Q r ee p r /(iM 2 y H 2 ), and 
fl* = p K /(3M 2 H 2 ) = -K/(aH) 2 . 



pi 

The autonomous equations for the variables ri, r 2 , Q r , and f2i<- are given by |98l 125 

r[ = ^(ri-l)ri[ri(ri(-3a + 4j8-2) + 6a-5j8)-5j9] 

x [2 (f2 r - Q K + 9) + 3r 2 (r?(-3a + 4/3 - 2) + 2r 2 (9a -90 + 2)- 15r x a + 14/3)] , (152) 
r 2 = -^-[r 2 (6r?(r 2 (45a 2 - 4(9a + 2)f3 + 36/3 2 ) - (Q r - n K - 7) (9a - 9/3 + 2)) + r?(-2(O r - Q K + 33) 



x 



(3a -A/3 + 2)- 3r 2 (-2(201a + 89)/3 + 15a(9a + 2) + 356/3 2 )) - 3na(-28fi r + 28n K + 123r 2 /3 + 36) 



+10/3(-lia- + 110 K + 21r 2 /3 - 3) + 3r 4 r 2 (9a 2 - 30a(4/3 + 1) + 2(2 - 9/3) 2 ) + 3r^r 2 (3a - 4/3 + 2) 2 

+3r 1 r 2 (9a - 9/3 + 2)(3a - 4/3 + 2))] , (153) 

Q' r = fl r (-A-2H'/H) , (154) 

Q' K = Q K (-2-2H'/H) , (155) 



5 In this case the only solution in the Minkowski background (H = 0) corresponds to <j> = for C2 ^ 0. 
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where a prime represents a derivative with respect to N = In a, and 

A ee 2r 4 r 2 [72a 2 + 30a(l - 5/?) + (2 - 9/3) 2 ] + 4r 2 [9r 2 (5a 2 + 9a/3 + (2 - 9^/3) + 2(9a -9/3 + 2)] 

+4r?[-3r 2 (-2(15a + l)/3 + 3a(9a + 2) + 4/3 2 ) - 3a + 4/3 - 2] - 24r ia (16r 2 /3 + 3) + 10/3(21r 2 /3 + 8) . (156) 

The Hubble parameter follows from the equation H' /H = — 5r[/(4ri) — r 2 /(4r 2 ). The solutions to Eqs. (| 154|) and 

(HHJ are given by Q r (JV) = fl^ e~ 4N H% / H 2 (N) and ^(iV) = Cl^ e~ 2N H% / H 2 {N) respectively, where the subscript 
"(0)" represents the values today (N = 0). 

From Eqs. (|152[) and (|153|) we find that there are three distinct fixed points: (A) (n, r 2 ) = (0, 0), (B) (n, r 2 ) = (1, 0), 
and (C) (ri,r 2 ) = (1,1). As we have already mentioned, the point (C) corresponds to the de Sitter solution. By 
considering homogeneous perturbations about this point, we can show that the de Sitter solution (C) is always 
classically stable [98[ . The point (B) is a tracker solution found in Ref. [§3], along which the field velocity evolves 
as (f> oc 1/H. During the radiation and matter eras the variable r 2 is much smaller than 1. The fixed point (B) is 
followed by the stable de Sitter point (C) once r 2 grows to the order of 1. If the initial conditions of both n and r 2 
in the radiation era are much smaller than 1, then the solutions are close to the point (A) at the initial stage. At 
late times the solutions approach the tracker at r% = 1. Depending on the initial values of n, the epoch at which 
the solutions reach the tracker is different. In the following we consider the background evolution in two regimes: (i) 
7"i = 1 and (ii) r\ <C 1 in more detail. 



1. Tracker solution (ri = 1) 

Along the tracker (n = 1) the dark energy density parameter (|151l) is given by 

Ode = r 2 , (157) 

which is much smaller than 1 during the radiation and matter eras. From Eqs. (|154|) and (|155|) we obtain r 2 /r 2 = 
8 + 2Q' r /Q r . This is integrated to give 

r 2 = d ia s n 2 r , (158) 
where d\ is a constant. From Eqs. ()153[) and (|154|) we have VI' k /Qk — £l' r /Clr = 2, which is integrated to give 

~f = d 2 a 2 , With d 2 = -Jjrr . (159) 

Substituting Eqs. (|158l) and f|159[) into Eq. (|154[) , we obtain the cosmologically viable solution to f2 r , as 



-1 + d 3 a - d 2 a 2 + V4dia 8 + (-1 + d 3 a ~ d 2 a 2 ) 2 

C2r = 771— g . ( 160 

where d 3 is another constant. Since the density parameter (I160p evolves as Q r ~ 1 + d 3 a in the early time (a <C 1), 
this demands the condition d 3 < (provided Ode > 0). Using the density parameters today, the constants d\ and 
d 3 can be expressed as d\ = [1 — — o[ 0) — f2^]/(r2r ') 2 and d 3 = — £1$ Since Cl r oc p r /H 2 oc l/(a 4 i? 2 ), 
we have H 2 /Hq — (£1^ /f2 r )(l/a 4 ). Using Eq. ()160|) . the Hubble parameter can be expressed in terms of the redshift 
z=l/a- 1: 

(^) 2 = ^(l + z) 2 + i^)(l + z)3 + I^)(l + ,) 4 

+] Ji n£> - o< 0) - nj°> + ii±^ [^i 0) + n£>(i + z) + n<°>(i + z)2] 2 , (lei) 

which is useful to test the viability of the tracker solution from observations. 
On the tracker, the equations of state defined in Eq. (|147|) are given by 

O r - CL K + 6 n r -Cl K - 6r 2 

»de = 57 — p— , Weff = 57 — rr • (162 

3(r 2 + 1) 3(r 2 + 1) 
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Figure 4: Evolution of uide versus z for a — 0.3 and /3 = 0.14 [cases (a)-(d)] in the fiat FLRW Universe (K — 0). We 
choose four different initial conditions: (a) n = 5.000 x 10~ n , r 2 = 8.000 x 10 -12 , and Q r = 0.999995 at z = 5.89 x 10 s , (b) 
n = 1.500 x 10" 10 , r 2 = 2.667 x 10" 12 , and Q r = 0.999992 at z = 3.63 x 10 8 , (c) n = 5.000 x 10" 9 , r 2 = 8.000 x 10" 14 , and 
Q r = 0.99995 at z = 6.72 x 10 7 , (d) n = 5.000 x 10" 6 , r 2 = 8.000 x 10" 17 , and ft r = 0.9986 at z = 2.04 x 10 6 . The case (e) 
corresponds to the tracker solution. From Ref. [Hj]. 



During the cosmological sequence of radiation (f2 r ~ 1, \Qk\ *C 1, r% <C 1), matter (Sl r <C 1, \^k\ <C 1, r 2 <C 1), and 
de Sitter (f2 r <C 1, \Qk\ <C 1, r 2 = 1) eras, the dark energy equation of state evolves as »de = —7/3 — > — 2 — > — 1, 
whereas the effective equation of state evolves as w c ff = 1/3— >0— ^— 1. This peculiar evolution of wde for the tracker 
corresponds to the case (e) in Fig. Although the effect of the cosmic curvature does not affect the dynamics of Wde 
significantly, it can change the diameter distance as well as the luminosity distance relative to the flat Universe. 

The epoch at which the solutions reach the tracking regime r% ~ 1 depends on model parameters and initial 
conditions. The approach to this regime occurs later for smaller initial values of r%, see Fig. 2J In Ref. [98[ it was 
shown that the tracker is stable in the direction of r\ by considering a homogeneous perturbation 6r%. This means 
that once the solutions reach the tracker the variable r% does not repel away from 1. If r\ < 2 initially, numerical 
simulations show that the solutions approach the tracker with the late-time cosmic acceleration. Meanwhile, for the 
initial conditions with r\ > 2, the dominant contribution to JIde comes from the Lagrangian C 2 , so that the field 
energy density decreases rapidly as in the standard massless scalar field. 



2. Solutions in the regime n <C 1 

There is another case in which the solutions start to evolve from the regime r\ <C 1 (where the term £5 gives the 
dominant contribution to the field dynamics) . In this regime, the variables n and r 2 satisfy the following approximate 
equations 

, 9 + n r - n K + 21/3r 2 , 3 + im, -H^-21/3r 2 

ri * sT^ ri ' r2 " sT^ r2 ■ (163) 

As long as {/3t" 2 , |^k|} *C 1, the evolution of T\ and r 2 during the radiation (matter) era is given by 7*1 oc a 5 ^ 4 and 
r 2 oc a 7 / 4 (n oc a 9 / 8 and r 2 oc a 3 / 8 ). Then the field velocity grows as </> oc t 3 / 8 during the radiation era and <j) oc i 1 / 4 
during the matter era. The evolution of <f> is slower than that for the tracker (i.e. <j> oc t). 
In the regime r\ <C 1 the equations of state are 

l + Q r -n K &n r - 80 K - 21/3r 2 .^^ 

^de — g + 21/3r2 , ^ ff ^ 3 (8 + 21/3r 2 ) ■ (164) 



26 



Provided that {/3r 2 , *C 1, one has wde — - 1/4, tt>eff — 1/3 during the radiation era and wde — — l/8,iu c ff — 

during the matter era. This evolution of u>de is quite different from that for the tracker solution. 

In Fig. 0] the variation of wde is plotted for a number of different initial conditions with r± <C 1 [which correspond 
to the cases (a)-(d)]. As expected, the solutions start to evolve from the value u>de — —1/4 in the radiation era. For 
larger initial values ofn they approach the tracker earlier. This tracking behavior also occurs in the presence of the 
cosmic curvature K 12591 . 



3. Conditions for the avoidance of ghosts and Laplacian instabilities 



Let us find a model parameter space in which the appearance of ghosts and instabilities can be avoided in covariant 
Galileon gravity. In doing so, we need to study a linear perturbation theory on the FLRW background. For simplicity 
we focus on the flat Universe with K = 0. Let us consider the perturbed metric 



ds 2 = -[1 + 2^(t,x)]dt 2 +d tX (t,x)dtdx l +a 2 (t)[l + 2$(t,x)}dx 2 



(165) 



where vp, $, and x are scalar metric perturbations. We have chose n the gauge 8cj> — without a non-diagonal scalar 
perturbation in the spatial part of the metric, i.e. d%jl = [l8l| . Taking into account two perfect fluids with the 
equations of state Wi — Pi/ Pi (i = 1,2), there are three propagating scalar degrees of freedom. The velocity potentials 
Vi (i — 1, 2) of perfect fluids are related with the energy-momentum tensor T^ l \ as Tj = —{pi + P^djVi (i = 1, 2). 
Introducing the vector Q = (vi,u 2 , < I>) and expanding the action (|140l) up to the second-order, we obtain the 



second-order action for scalar perturbations [H[ (see also Refs. [l83l l26Cll . l26ll |): 



SS 



(2) 



1 



dtd 3 xa 3 



Q*AQ- 



1 



rVQ'CVQ- Q'BQ- Q l DQ 



(166) 



where the fields "J and x are integrated out. A, C and D are 3x3 symmetric matrices and B is an antisymmetric 
matrix (for which we do not write explicit forms). 

In order to avoid the appearance of ghosts we require that the matrix A is positive definite. This corresponds to 
the conditions (1 + wi)pi/wi > 0, (1 + W2)pijw2 > 0, and 



Q s _ 6(1 + pi)(pi + p 2 + M1M2 - 2/i 3 - pi) 



(i+p 3 y 



> 0, 



where 



pi = 3arir 2 /2 - 3/3r 2 , 

p 2 = (3a - 4/3 + 2)r?r 2 /2 - 2(9a -9/3 + 2)r 2 r 2 + 45ar!r 2 /2 - 28/3r 2 , 
p 3 = -(9a - 9/3 + 2)r 2 r 2 /2 + 15arir 2 /2 - 21/3r 2 /2 . 

The propagation speeds cs of three scalar degrees of freedom are known by solving the equation 

det(c|A - C) = . 



(167) 



(168) 
(169) 
(170) 



(171) 



For the two perfect fluids we have c 2 s = w\ and c 2 s = w 2 , which are are positive for both radiation and non-relativistic 
matter. The third stability condition associated with another scalar degree of freedom is given by 



C S 



where 



(1 + Mi) 2 [2/4 - (1 + M3X5 + 3w cff ) + m r + 3Q m ] - 4^(1 + + p 2 ) + 2(1 + p 3 ) 2 (l + p A ) 

6(1 + Pl){pi + P2 + P1P2 - 2p3 - p\) 



>0, (172) 



Pi = -arir 2 /2- 3^r 2 (ri/ri +r 2 /r 2 )/4. 



Hi = a 2 hij, where hij is traceless (h l 



Let us consider tensor perturbations with 
0). We expand the action (|140p at second-order in terms of the two polarization modes, h tJ 
e® and are the polarization tensors. For the polarization mode /i®, the second-order action is given by 



(173) 

0) and divergence- free j — 

i @ e® + ftg, e|, where 



dt d 3 x a 3 Q T 



(174) 



27 



a= 12 i 3 l,2 -9 l 8-2, fr 2 = 0at(c T 2 )' 



a = 2/1+ j, (&), i=1 =0 




a = 2ft(c/) ds =0 
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13 13 ls »' 





Figure 5: The viable parameter region (blue colored region) along the tracker solution (n = 1) in the covariant Galileon model 
(where r-z > 0). If the solutions start from the regime r\ <C 1, we also require that /3 > 0. From Ref. [971 ]. 



The conditions for the avoidance of ghosts and Laplacian instabilities of tensor perturbations correspond, respectively, 
to 



= - + -arir 2 - 7r/3r 2 > . 



Qt 

Ml . , 
c 2 = 2ri(2 - ar 1 r 2 ) - 3f3(r 2 r[ + rir' 2 ) Q 
Ct ~ 2ri(2 + 3arir 2 -6/3r 2 ) 



(175) 
(176) 



The same conditions also follow from h®. 

In the regime r\ -C 1 and r 2 -C 1 one has Qg/M^ ~ 60/3r 2 and — 1/2- For the initial conditions with r 2 > 
we require that (3 > to avoid the scalar ghost. Since c| ~ (1 + O r )/40 and cf, ~ 1 + 3/3r 2 (5 — 3fi r )/8 ~ 1, there are 
no Laplacian instabilities of scalar and tensor perturbations in this regime. 

In the tracking regime characterized by r% = 1 (either r 2 <C 1 or r 2 = 1), the conditions (I167p . (I172p . (I175p . and 
(11761) give the bounds on the parameters a and j3. In the regime r 2 <C 1 these conditions translate to 



4 



s ~ 3(2 - 3a + 6/3)r 2 > , 

! + 10a - 9£ + n r (2 + 3a - 3/3) 



3(2 -3a + 6/3) 



> 0. 



(177) 
(178) 



For the branch r 2 > the first condition reduces to 2 — 3a + 6/3 > 0. Since c\ ~ 1 — r 2 (4a + 3/3 + 3/30 r ) /2 ~ 1 and 
Qr/Mpj = 1/2 + 3(a — 2/3)r 2 /4 > 0, the tensor modes do not provide additional constraints. At the de Sitter point 
(r 2 = 1) we require that 



Qs 4-9(a-2/3) 2 



3(a - 2/3) 2 



> 0, 



Qt 
Ml 



1 



2 _ (a - 2/3)(4 + 15a 2 - 48a/3 + 36/3 2 ) 
Co — — t, —, t-^t-, > . 



(2 + 3a - 6/3) > , 



2[4-9(a-2/3) 2 



2 + 3a - 6/3 



> 0. 



(179) 
(180) 



If /3 > 0, it can happen that c?n has a minimum during the transition from the regime r 2 -C 1 to r 2 ~ 1 (97l . [96 
This value tends to decrease as /3 approaches 1. Imposing that > at the minimum, we obtain the bound 



a < 12^/3-9^-2. 



(181) 



In Fig. [S]we plot the parameter space in the (a,/3) plane constrained by the conditions (|177[) - (|181[) . Clearly there are 
viable model parameters satisfying all the theoretical constraints. 
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4- Observational constraints on Galileon cosmology from the background cosmic expansion history 

Since the evolution of the dark energy equation of state in covariant Galileon gravity is rather peculiar, the ob- 
servational data related with the background cosmic expansion history may place tigh t con straints on the model. 
Especially the analytic formula fllOlf) for the tracker is useful for such a purp ose. I n Ref. |259j the au thor s confronted 
the Galileon model by using the observational data of SN la (Constitution |262| and Union2 sets |263| ). the CMB 
(WMAP7) shift parameters Q, and B AO (SDSS7) @. 

If either of the SN la data (Constitution or Union2) is used in the data analysis, the % 2 for the tracker is similar 
to that in the ACDM model. In the presence of the cosmic curvature K, the tracker solution is compatible with 
the individual observational bound constrained from either CMB or BAO. However, the combined data analysis of 
Constitution+BAO+CMB shows that the difference of \ 2 between the tracker and the ACDM is 5\ 2 ~ 22 (or ~ 4.3a). 
This means that the tracker is severely disfavored with respect to the ACDM. A similar conclusion was reached from 
the combined data analysis of Union2+BAO+CMB. The reason for this incompatibility is that the SN la data favor 
the large values of ®f£) (^ 0.32), whereas the CMB and BAO data constrain smaller values of 0$ 0.27). 

The general solutions starting from the regime r\ -C 1 finally approach the tracker as T\ grows to 1. In Ref. |259l ] 
the authors carried out the likelihood analysis for such general solutions and found that the solutions approaching the 
tracker at late times (such as the case (a) in Fig. [4]) are favored from the combined data analysis. In the flat FLRW 
background the best-fit model parameters are a = 1.411 ±0.056, (3 = 0.422 ±0.022 (Constitution±CMB+BAO, 68% 
CL), and a = 1.404 ± 0.057, (3 = 0.419 ± 0.023 (Union2+CMB+BAO, 68% CL). 

For several fixed values of it was shown that the late-time tracking solutions can be consistent with the data, 

apart from the models with largely negative fl^P such as < —0.01. For example, the general solutions with 

= 0.01 and the model parameters {a, (3) = (1.862,0.607) g ive t he similar value of \ 2 to that in the the ACDM. 
In this case t he A kaike-Information-Criteiron (AIC) statistics |264j | also have the same support for the two models 
(see also Ref. [265] ) 

The Bayesian-Information-Criterion (BIC) statistics [266j | show that the general solutions, with all 4 parameters 

(a, /?, Q,m\ ^jP) are varied, are not particularly favored over the ACDM model. This mainly comes from the statistical 
property that the numbers of model parameters are larger than those in the flat ACDM. In fact the late-time tracking 
solutions with a non-zero cosmic curvature can be well consistent with the combined data analysis at the background 
level. 



B. Generalized Galileon gravity 

In Sec. IIIII we showed that in Brans-Dicke theory with the coupling Q of the order of unity the presence of the field 
potential allows a possibility for the consistency with local gravity constraints through the chameleon mechanism. 
Another way to recover the General Rel ativis tic behavior in the regions of high density is to introduce the Galileon-likc 
field self-interaction. Silva and Koyama [238| studied Brans-Dicke theory in the presence of the term ^((/))D(j)(d '^tfid' 1 (f>) 
[which is the generalization of the term £3 = O^d^tfid^ 4>)]. The action of this theory is given by 



S = 



(V0) 2 + £(0)Q^<W) 



d 4 xC 



M ■ 



(182) 



If £(0) oc <fi~ 2 , there exists a de Sitter solution that can be responsible for the late-time acceleration. As in the Galileon 
model discussed in Sec. IV Al the field is nearly frozen during the radiation and matter eras through the cosmological 
Vainshtein mechanism, but it finally approaches the de Sitter solution characterized by <\> = constant. Moreover, as in 
the DGP model, the Vainshtein radius can be much larger than the solar system scale, so that the General Relativistic 
behavior can be recovered in the local region |238| . 

We may consider more general theories described by the action |242| 



S = 



-F{^)R 



B(<f>)X + ^U^d^d^) 



d 4 xC 



M ■ 



(183) 



where X = —g fJlI/ d fl (j>d u cf)/2, and F((j>), B(<f>), £(</>) are functions of 4>. From the requirement of having de Sitter 
solutions responsible for dark energy, it is possible to restrict the functional forms of F(<fi), B((f>), and £(</>). In the 
presence of non-relativistic matter (energy density p m ) and radiation (energy density p r ), the field equations are given 
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by 

6F F + F { 6£ + 3FH 2 + 3FH 2 1 



_ 2 H_ = (*F± _ 2^ 
H 2 \ F F 




Btf F^ 2 <j>F A 6^ 3 2 2 ^ 

■x H — x — I — — H x 1 



F F F F \ 3£ 



FH 2 3FH 2 



(185) 



where x = <j>/(H<p). Let us search for de Sitter solutions at which iJ and x are constants. If F(0) and £(</>) are 
power-law functions of <j>, the quantities such as cfrF^/F, F^</> 2 /F, and remain constants. Provided that 

B/F oc 0~ 2 and £/F cx <fi~ 3 , we can solve Eqs. (|184[) and (|185|) for a; and -ff at the de Sitter point. These conditions 
are satisfied for the following functions 

F(4>) = MpV^/Mp!) 3 "" , B(<t>) = w^/MpO 1 "" , ^(0) = (A//z 3 )(4>/M pl )-" , (186) 

where M p \ ~ 10 18 GeV is the reduced Planck mass, /i (> 0) is a constant having a dimension of mass, and ui and A are 
dimensionless constants. One can show that the coupling A must be positive for the consistency of theories [242] . The 
Brans-Dicke theory described by the action (|182[) corresponds to n = 2 with the Brans-Dicke parameter wbd = 
Since F(<f>) is constant for n = 3, the theory with n — 3 corresponds to k-essence minimally coupled gravity. 
From Eqs. (|184[) and (|185[) we obtain the following algebraic equations at the dc Sitter fixed point: 



n(n - 3) 2 x 3 dS + (n - 3)(n - 12)x^ s - 6(n - 5):r d s + 18 

p 3 [(n - 3)x ds - 2}[(n - 3)x dS - 3} 
M pl H 2 s 2x 3 s (x dS +3) 



u = ' ^ ^ f-f i ' ™ ' , (187) 

^ds^s + 3) 

(188) 



where XdS an d -HdS are the values of x and H at the de Sitter point, respectively. We fix the mass scale p to be 
p = (MpiiJjg) 1 / 3 , where we have used -ff d s — 10 _60 M p i. For given w and n, the quantity x^s is determined by solving 
Eq. (|187j) . Then the dimensionless constant A is known from Eq. (|188[) . 

In order to recover the General Relativistic behavior in the early cosmological epoch we require that the field initial 
value 4>i is close to M p \ from Eq. (II 86[) . The quantity x is much smaller than 1 in the early cosmological epoch, so that 
the field is nearly frozen during the radiation and matter eras. The field starts to evolve at the late cosmological epoch 
in which x grows to the order of unity. Introducing the dimensionless quantities y = Xx 2 H 2 / H dS and Q r = p r I (3 FH 2 ) , 
one can show that the fixed point corresponding to the matter era corresponds to (x, y, Q r ) = (0, (3 — n)/6, 0) j242j |. 
Since y is positive definite, it follows that n < 3. 

The conditions for the avoidance of ghosts and Laplacian instabilities are known by employing the method presented 
in Sec. IV A 31 Provided F(0) > 0, we require that x > to avoid ghosts during the cosmological evolution from the 



radiation era to the epoch of cosmic acceleration 242] . The stability of the de Sitter point is automatically ensured 



for n < 3 and x^s > 0. From Eq. (|187[) the parameter u is restricted in the range 

ui < -n(n - 3) 2 . (189) 

The field propagation speed squared during the radiation and matter eras is given by c| ~ 6/5 and c 2 s ~ 2/3, 
respectiv ely i n which case no instabilities of linear perturbations are present. Meanwhile, at the de Sitter solution, 
we have [242j ] 

c 2 = (» - 2)[(n - 3)(n - 4)x 2 s - 8(n - 3)z ds + 6]^ dS , , 

(n - 3) (3n 2 - lOra + 12) a;^ + 12(n - 5)x 2 s + 18(n - 8)a;ds - 108 ' { 1 

which is positive for n > 2. Hence the parameter n is restricted in the range 

2<n<3, (191) 

which includes Brans-Dicke theory with the action (|182[) as a specific case (n — 2). 

In Ref. [239j the authors studied the evolution of matter density perturbations and showed that, for the model with 
n = 2, there is an anti-correlation between the cross-correlation of large scale structure and the integrated Sachs- Wolfe 
effect in CMB anisotropies. We shall discuss the main reason of this anti-correlation in Sec. IVII Dl This property will 
be useful to distinguish the above model from the ACDM in future observations. 
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VI. OTHER MODIFIED GRAVITY MODELS OF DARK ENERGY 

In this section we briefly discuss other classes of modified gravity models of dark energy. These include (i) Gauss- 
Bonnet gravity with a scalar coupling f{<j>)G, (ii) R/2 + f(Q) gravity, and (iii) Lorentz- violating models. 

A. Gauss-Bonnet gravity with a scalar coupling 



In addition to the Ricci scalar R, we can construct other scalar qu antit ies coming from the Ricci tensor and 
the Riemann tensor R^ap, i.e. P = R tlu R ,JjV and Q = R,,„r,nR^ va ^ It is po ssible to avoid the appearance of 

spurious spin-2 ghosts by taking a Gauss-Bonnet (GB) combination |268t I269L l27l| . defined by 



G = R 2 - iR^RP" + R^apR!""* . 
A simple model that can be responsible for the cosmic acceleration today is (4E 



d 4 2H 



1 



:R 



1 



(v^) 2 - v(<j>) - fi^g 



d 4 xC 



M 



(192) 



(193) 



The coupling of the field with the GB term appears in low- 



For the exponent ial potent ial V{<f>) = Voe A * and the coupling f(<f>) — (f //j,)e 



where V(4>) and f(4>) are fun ction s of a scalar field 
energy effective string theory [270 

the cosmological dynamics were studie d in Refs . (43l. [282h428 7j . In this model there exists a de Sitter solution due to 
the presence of the GB term. In Refs. [282l |283| | it was also found that the late-time de-Sitter solution is preceded by 
a scaling matter era. 



Koivisto and Mota 282] placed observational constraints on the model (|193[) with the exponential potential V((j>) = 
Vbe~ A< ^, by using the Gold data set of SN la together with the CMB shift p aram eter data of WMAP. The parameter 
A is constrained to be 3.5 < A < 4.5 at the 95% confidence level. In Ref. 284], they also included the constraints 
coming from the BAO, LSS, big bang nucleosynth esis, and solar system data. This joint analysis showed that the 
model is strongly disfavored by the data. In Ref. |283j | it was also found that tensor pe rturbations are subject to 
negative instabiliti es wh en the GB term dominates the dynamics (see also Refs. [288l |289| for related works). 

Amendola et al. 290] studied local gravity constraints on the above model and showed that the density parameter 
r^GB coming from the GB term is required to be strongly suppressed for the compatibility with solar-system experi- 
ments (which is typically of the order of J7gb < 10~ 30 ). The above discussion indicates that the GB term with the 
scalar-field coupling f(4>)G can hardly be the source for dark energy. 



B. R/2 + f(Q) gravity 



The general Lagrangian including the scalar quantities constructed from the Ricci scalar R, the Ricci tensor R^ v 
and the Rieman n tensor R {Lvn R is given by C = f(R, P, Q). The dark energy models based on these theories have been 
studied in Refs. I272ll279l. In order to avoid spurious spin-2 ghosts we need to choose the GB combination (|192p . i.e. 
C = f(R, Q - 4P) [280] [281|. 

The cosmological dynamics based on the action, 



S 



<FxC 



M 



(194) 



have been studied by a number of authors [44ll29ll4296j ]. In order to ensure the stability of radiation/matter solutions, 
we n eed t o satisfy the condition f : gg > for all g. We also require the regularities of the functions /, f t g, and 
f gg 293]. There exists a de Sitter solution responsible for dark energy, whose stability requires the condition 
< HdS.f. Ga(H ds) < 1/384. Moreover f.gg must approach +0 in the limit \g\ —> oo. 

In Ref. [293] the authors proposed a number of /((?) models satisfying these conditions (see also Ref. [294j ). One 
of such models is given by 



f(G) = A- 



, i G 

, zarctan — 

rgi \g* 



(195) 



where a, A and G* are constants. The numerical simulation of Ref. 293] shows that the model (I195[) is cosmologically 
viable at least at the back ground level. Moreover it can be consistent with solar-system constraints for a wide range 
of model parameters |297l ] . 
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In order to study the viability of the theories described by the action (|194|) further, let us consider the evolution of 
matter density perturbations in the presence of a perfect fluid with the barotropic equation of state % = Pm/pm- 
In Ref. [298| it was shown that, for small scales (i.e. for large momenta k), there are two different scalar propagation 
speeds. One of them corresponds to the mode for the fluid, i.e. c 2 = wm, whereas another is given by 

2 2H l + w M K 2 p M , . 

C2 = 1+ iP + TW^ F ' (96) 



where 



M = HQfm ■ (197) 



The parameter fi characterizes the deviation from the ACDM model (note that the linear term f = c Q do es not give 
rise to any contribution to the field equation). For viable f(G) models we have <gc 1 at high redshifts |293j . Since the 
background evolution during the radiation/matter domination is given by 3H 2 ~ k 2 Pm and H / H 2 ~ — (3/2)(1 + u>m), 
it follows that 

c£ ~ -1 - 2w M ■ (198) 

The Laplacian instability at small scales is absent only for wm < — 1/2. Since wm = 1/3 and wm = during the 
radiation and matter eras, respectively, the perturbations with large momentum modes are unstable. Th is leads t o 
violent growth of matter density perturbations incompatible with the observations of large-scale structure |292l |298| | . 

By co nsidering the full perturbation equations, one can show that the onset of the negative instability corresponds 
to (298| 

H « {aH/kf . (199) 

Even when fi is much smaller than 1, we can always find a wave number k (3> aH) satisfying the condition (|199l) . 
For the scales smaller than that determined by the wave number in Eq. (11991) . the linear perturbation theory breaks 
down. Hence the background solutions cannot be trusted for those scales, which makes the theory unpredictable. The 
Laplacian instability can be avoided only for fi — 0, which corresponds to the ACDM model. The above property 
persists irrespective of the forms of f{Q). 

For more gene ral theories described by the Lagrangian dens ity f(R, G) it is possible to avoid such Laplacian 
instabilities [299j |. depending on the models [300] (see also Refs. |30ll4304| ). It may be of interest to construct some 
viable dark energy models in such theories. 



C. Lorentz violating models 

The modified gravity models such as f(R) gravity and Galileon gravity can give rise to the phantom dark energy 
equation of state wde < — 1 without violating the conditions for the appearance of ghosts and instabilities. In the 
models with a broken Lor entz invariance it is also poss ible to realize wde < — 1 without pathological behavior in the 
Ultra- Violet (UV) region [305M307| (see also Ref. 308] for a review). In order to construct Lorentz violating models 



without pathological behavior of phantoms, one may start with a field theory consistent at energy scales from zero 
to the UV cutoff scale Ai and then deform the theory in the Infra-Red (IR) in such a way that its behavior at high 
energies remains healthy. Although the weak energy condition is violated in the homogenous background, pathological 
states are present below a certain low scale e only. Provided that e is close to the Hubble scale, a theory of this sort 
should be acceptable. 

L et us consider a Lorentz violating model with two-derivative kinetic term s with healthy behavior below the scale 
Ai [3091 ] plus one-derivative term suppressed by the small parameter e |310| . The model has a vector field and a 
scalar field $ with a potential V(B, $). The Lagrangian is given by [306l l307j 



C = C {2) + + £ (0) + C M , (200) 
C {2) = -^a{~)g» x D^B„D^B x + ^(SJ^B^Ba^- + \d^d^ , (201) 
£« = ed^B" , C {0) = -V(B, $) , (202) 



where S = B^B^ /Ai 2 (Ai is the UV cut-off scale). The dimensionless parameters a and (3 are the functions of S, 
and e is a free positive parameter that characterizes an IR scale. The Lorentz invariance is broken for H ^ 0. 
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Figure 6: Evolution of u>de (denoted as w in the figures) and w e s = 1 — 2H /(3H 2 ) together with the density parameters 
of dark energy (Qde), non-relativistic matter (f2 m ), and radiation (f2 r ) for the Lorentz-violating model (|200|) with V = 
m 2 </> 2 /2 — M 2 X 2 jl. The model parameters are chosen to be a = 1, 7 = 1/2, e/m = 3 and M/m — 1 with the initial conditions 
X = X = <f> — 0, 4> — 0.50/1, and f2 r = 0.99, TO = 0.01 [where (f>A = Mm p i / (V 4ftam) ] . The present epoch corresponds to 
fi m = 0.3 and Qde = 0.7, which is denoted by a vertical line. After the cosmological constant boundary c rossin g, idde reaches 
a minimum —1.19 and then it increases toward the de Sitter value —1 from the phantom side. From Ref. 307]. 



In spatially homogeneous background we have 

B = X , 



Bi = , 



(203) 



where Bq and Bi are time and space components of B^. In the flat FLRW background with non-relativistic matter 
and radiation for the matter Lagrangian C m , we obtain the following equations of motion 



H- 



8ttG 



7 [X 



3 

3HX 



l 7 X 2 - ^H 2 X 2 + i^ 2 + W(4>) + U(X) + Pm + p r 



1 



7, x i 2 + -a :X H 2 X 2 + ZaH 2 X - ej> + V x = , 



3H<j) + e(X + 3HX) + V,^ = , 



(204) 

(205) 
(206) 



where j(X) = X 2 f3(X)/M 2 - a(X). 

For the separable potential V — m 2 cj) 2 1 2 — M 2 X 2 /2 with constant X the cosmological dynamics of the above system 
have been studied in detail in Refs. |306l 1307], In what follows we assume that both a and 7 are constants. Provided 
e/m > y // 2a/3 one can show that there is a de Sitter solution with H = M/y/Sa, at which and X are frozen. In the 
early cosmological epoch the field X is close to 0, whereas the field (f> slowly rolls down its potential (with the energy 
density dominating over that of X). 

Prior to the epoch of de Sitter cosmic acceleration there is a transient phantom regime (iude < — 1) characterized by 
H < Mj \fioL in which the field <f> rolls up its potential. The Hubble parameter slowly increases toward the de Sitter 
value if ds = M/y/3a. In Fig. [6] we plot the evolution of the dark energy equation of state u>de = -Pde/pde and the 
effective equation of state w c s = 1 — 2H/(3H 2 ) together with the density parameters of dark energy, non-relativistic 
matter, and radiation 6 . Clearly the solution undergoes the period with u>de < —1 before reaching the de Sitter 
attractor. 



Here the definition of p DE and P DE is p DE = yX 2 /2 + <j) 2 /2 + V - 3aH 2 X 2 /2 and P DE = yX 2 /2 + <j) 2 /2 - V + eif>X + aHX 2 + 
2aHXX + 3aH 2 X 2 /2. 
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The phantom equation of state can be realized without having ghosts, tachyons or superluminal modes in the 
UV region. In the IR region characterized by p < e, either tachyons or ghosts appear for the spatial momenta p 
smaller than y/ (e 2 — M 2 )/a — m 2 306, 307]. The presence of tachyons at IR scales leads to the amplification of 
large-scale field perturbations whose wavelengths are roughly comparable to the present Hubble radius. There are 
two tachyonic regions of spatial momenta in this model: (a) one is sub-horizon and its momenta are characterized by 
M 2 /7 < p 2 < (e 2 — M 2 )/a — m 2 ; (b) another is super-horizon with the momenta < p 2 < m 2 M 2 /e 2 . In the region 
(a) there is a parameter space in which the perturbations always remain smaller than the homogenous fields. In the 
region (b) the growth of the perturbations is suppressed by the factor m 2 /e 2 . 

There are ot her c l asses of Lorentz violating models such as ghost condensate [305j ] and Horava-Lifshitz gravity [31 lj ] 
(see also Refs. l312l| - |317| ). The application of such scenarios to dark energy has been studied by a number of authors, 
see e.g., [318l - l320j . It remains to to see whether such Lorentz- violating models can be observationally distinguished 
from other dark energy models. 



VII. OBSERVATIONAL SIGNATURES OF MODIFIED GRAVITY 



In order to confront modified gravity models with the observations of large-scale structure and CMB, we discuss 
the evolution of density perturbations in four modified gravity models: (i) f(R) gravity, (ii) scalar-tensor gravity, 
(iii) DGP braneworld model, and (iv) Galileon gravity. We also discuss observables to confront with weak lensing 
observations. 



A. f{R) gravity 



Let us first consider metric f(R) gravity in the presence of non-relativistic matter. We take the following perturbed 
metric in a longitudinal gauge about the fiat FLRW background with scalar metric perturbations <& and ^> [l8ll | 



ds 2 



(1 + 2*)di 2 + a 2 (l + 2<f>)6 ij dx i dx j 



(207) 



The energy momentum tensors of a non-relativistic perfect fluid are decomposed into background and perturbed parts, 
as Tq = ~(p m + Sp m ) and T° = —p m v m ,a (v m is a velocity potential). 



The equations for matter perturbations, in the Fourier space, are given by [321 323] 

5p m + 3H5p m = -p rn 3$ + {k 2 /a)v m , (208) 

v m + Hv m = tf/a , (209) 
where k is a comoving wave number. We define the gauge- invariant matter density perturbation 5 m , as 

S m = 8 p rn / p rn + 3Hv , where v = av m . (210) 



Then Eqs. (gOg) and (gD5| yield 



from which we obtain 



S m = -(k 2 /a 2 )v - 3($ - Hv)' 
v = * , 



2HS m + (fc 2 /o 2 )* = 3B + 6HB , 



(211) 
(212) 



(213) 



where B = — $ + Hv. 

In f(R) gravity the quantity F(R) — df /OR has a perturbation SF. In the following we use the unit k 2 — 8ttG = 1, 
but we restore gravitational constant G when it i s required. For the action given in Eq. @, we obtain the linearized 
perturbation equations in Fourier space |324| - |326| 



k 



1 



3H5F - ( 3H + 3H 2 - ^ ) 5F - 3HF^ - 3F(HV - $) - 5p n 



+ 3H(H9 -*) = _ 

fk 2 \ 1 • • 

SF + 3HSF + ( -J + M 2 J SF = -Sp m + F(3H^ + * - 3$) + (2F + 3HF)$> , 

* + $ = -5F/F . 



(214) 

(215) 
(216) 
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In Eq. (|215[) we have introduced the mass term 



M 2 = - [ -f R ) . (217) 



3 V/, 



,RR 



For viable dark e nergy models the condition Fj f t RR 3> R is satisfied during most of the cosmological epoch, so that 



M 2 ~ F/(3f tRR ) [69|,l7l|. This is equivalent to the mass squared M 2 introduced in Eq. (l4"3")l , which is required to be 



positive to avoid the tachyonic instability. 

For the observations of large-scale structure and weak lensing we are interested in the modes deep inside the Hubble 
radius (k 3> aH). In the following we employ the quasi-static approximation under which the d ominant terms in 
Eqs. (j2T5j) - ([2"il)l) correspond to those including k 2 /a 2 , 6p m (or 8 rn ) and M 2 [T^. [l55l . [l65l [327l - [329j . We then obtain 
the following approximate relations from Eqs. (|213|) - (j216[) : 

8 m + 2H8 m + (k 2 /a 2 )^ = 0, (218) 

* = 2b {P** - 5F )> * = -5F {P p ™ + SF ) • ^ 

8F + 3H8F + (k 2 /a 2 + M 2 ) SF = Sp m /3 . (220) 

The evolution of perturbations is different depending on whether M 2 is larger or smaller than k 2 /a 2 . We shall discuss 
two cases: (A) M 2 >• k 2 /a 2 and (B) M 2 <C k 2 /a 2 , separately. For viable f(R) models the mass squared M 2 is large 
in the past and it gradually decreases with time. Hence the transition from the regime (A) to the regime (B) can 
occur in the past, depending on the wave numbers k. 



1. Evolution of perturbations in the regime: M 3> k /a 

The solutions to Eq. (|220l) are given by the sum of the oscillating solution ^F osc obtained by setting 8p m = and 
the special solution 8F m< \ of Eq. (1220[) induced by the presence of matter perturbations 8p m . The oscillating part 
SF osc satisfies the equation (a 3 / 2 <5F osc ) " +M 2 (a 3,/2 (5Fosc) — 0. Using the WKB approximation, we obtain the following 
solution fsoj ) 

<5F osc ex a' 3 / 2 f. RR 1/4 cos ( / —2=dt) , (221) 

\J yij.RR ) 

where we have used the approximation F ~ 1. 

For the analytic estimation of the oscillating mode we take the model (|4"7) , which corresponds to the asymptotic 
form of the models (|45[) and (I46[) in the region R^$> R c . During the matter era in which the background Ricci scalar 
evolves as R (0) = 4/(3£ 2 ), the quantity f, RR has a dependence f, RR oc J R" 2 ( n+1 ) cx t 4(n+1 ) for the model (J47J). Then 
the evolution of the perturbation, SR OS c — 8F osc / f, RR , is given by 

<5i? osc ~ c<- (3 ™ +4) cos(c t- 2(rl+1) ) , (222) 

where c and Co are constants. Unless the coefficient c is chosen to be very small, as we go back to the past, the 
perturbation SR osc dominates over the background value i?(°)(oc t~ 2 ). Since the Ricci scalar can be negative, this 
leads to the violation of the stability conditions (f RR > and F > 0). 

The special solution SFi n d to Eq. (|220p can be derived by neglecting the first and second terms relative to others, 
giving 

<5F ind ~ S Pm /(3M 2 ) , SR ind ~ 8p m . (223) 

Under the condition |<5F osc | < |<5F ind |, one has 8F ~ 8p m /(3M 2 ) and hence * = -$ = -(a 2 /k 2 )8p m /{2F). Then the 
matter perturbation equation (|218|) reduces to 

8 m + 2H8 m - 4iTGp m 8 m /F = , (224) 



where we have reproduced the gravitational constant G. In Refs. [3301 133 lj the perturbation equations have been 
derived without neglecting the oscillating mode. 
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During the matter-dominated epoch (f2 OT = p m /(3FH 2 ) — 1), Eq. (|224l) has the growing- mode solution 

S m oc t 2 / 3 . (225) 

This is the same evolution as that in standard General Relativity. From Eq. (|223l) the matter-induced mode evolves 
as <5Fj n d oc t 4 ("+ 2 / 3 ) and <5i?i n d oc £~ 4 / 3 . Compared to the oscillating mode (j222[) . <5i?i n d decreases more slowly and 
hence it dominates over SR osc at late times. The evolution of the perturbation 5R = SR osc + 5i?ind relative to the 
background value R(°> is given by 

5R/RW ~ ci t-( 3 "+ 2 ) cos(cot^) + c 2 i 2/3 , (226) 

where c\ and c 2 are constants. In order to avoid the dominance of the o scilla ting mode at the early cosmological 
epoch, we require that the coefficient c\ is suppressed relative to C2 [69l. ItD . l332j | . 

This fine-tuning of initial conditions is related with the singularity problem raised by Frolov [333$ . The field 
<f> = ^/3/2k 2 \s\F in the Einstein frame has a weak singularity at <p = (at which the curvature R and the mass 
M go to infinity with a finite potential V). Unless the oscillating mode of the field perturbation 8<j> is strongly 
suppressed relative to the background field (f>^ , the syste m can access the curvature singularity. This past singu larity 
can be cured by taking into account the R 2 term [334j . Note that the f(R) models proposed in Ref. l335ll [e.g.. 
f{R) = R — aR c ln(l + R/R c )] to cure the singularity problem satisfy neither the local gravity constraints [336| nor 
observational constraints of large-scale structure 337J. There ar e some works for the construction of unified models 
of inflation and dark energy based on f(R) theories |34 1338ll339j | . but the smooth transition between two accelerated 
epochs without crossing the point J,rr — is not easy unless the forms of f(R) are carefully chosen |334| . 



2. Evolution of perturbations in the regime: M 2 <C k 2 /a 2 

Since the mass M decreases as M ex t~ 2tyn+1 \ the modes initially in the region M 2 3> k 2 /a 2 can enter the regime 
M 2 <C k 2 /a 2 during the matter-dominated epoch. It is sufficient to consider the matter-induced mode because the 
oscillating mode is already suppressed during the evolution in the regime M 2 3> k 2 /a 2 . The matter-induced special 
solution of Eq. (|220l) in the regime M 2 <C k 2 /a 2 is approximately given by 

a 2 

SF md ~ — ^ m . (227) 

From Eq. (|219[) the gravitational potentials satisfy 

4 1 a 2 r 2 1 a 2 ,„ . 

* = -3"2F^' *=3-2FF^- (228) 
Plugging Eq. (|228p into Eq. (|2 18[) , the matter perturbation obeys the following equation 

S m + 2H6 m - | • 47rGp m S m /F = . (229) 

During the matter-dominated epoch (f2 m ~ I and a oc i 2 ^ 3 ), we obtain the following evolution 

S m oc t (V33-i)/a . (230) 
The growth rate of 5 m gets larger compared to that in the regime M 2 3> k 2 /a 2 . 



3. Matter power spectra and the ISW effect 

If the transition from the regime M 2 k 2 /a 2 to the regime M 2 <C k 2 /a 2 occurs during the matter era, the 
evolution of matter perturbations changes from S m oc t 2 / 3 to S m oc ^(\ / 33-i)/6 use ^ e SUDSCr ipt "fc" for the 

quantities at which k is equal to aM, whereas the subscript "A" is used at which the accelerated expansion starts 
(a = 0). While the redshift z\ is independent of fc, Zk depend on k and also on the mass M. 

For the model (|47p the variable m = Rj \rr) f r can grow fast from the regime m (aH/k) 2 (i.e., M 2 ^> k 2 /a 2 ) to 
the regime m ^> (aH/k) 2 (i.e., M 2 <C k 2 /a 2 ). In fact, m can grow to the order of 0.1 even if m is much smaller than 
10 -6 in the deep matter era. For the sub-horizon modes relevant to the galaxy power spectrum, the transition at 
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Figure 7: The matter power spectra Ps m (k) in Brans-Dicke theory with the potential ((98}. This analysis covers the f(R) model 
(|47[l with the correspondence Q = — l/v6 and p = 2n/(2n + 1). Each case corresponds to (a) Q = 0.7, p = 0.6, C = 0.9, (b) 
Q = -1/VE, p = 0.6, C = 0.9, (c) the ACDM model, and (d) the ACDM model with a nonlinear halo-fitting (cr 8 = 0.78 and 
shape parameter V = 0.2). The model parameters are = 0.28, Ho = 3.34 x 10~ 4 /iMpc" 1 , n a = 1 and = 3.2 x 10~ 10 . 
From Ref. [34Pj ] . 



M 2 = A: 2 /a 2 typically occurs at the redshift £/. larger than 1 (provided that n — 0(1))- For the mode k/ (aoH ) = 300 
one has Zk = 4.83 for n = 1 and Zk = 2.49 for n = 2. As n gets larger, the period of non-standard evolution of S m 
becomes shorter because Zk tends to be smaller. Since the scalaron mass evolves as M oc i _2 (™+ 1 ) for the model (|47|) . 
the time tk has a scale-dependence tk oc k~ 3 /( 6n + 4 ) . This means that the smaller-scale modes cross the transition 
point earlier. The matter power spectrum Pg m — \S m \ 2 at the time t\ shows a difference compared to the case of the 



ACDM model: 



Ps m (t A ) ft A y\ b °J 



Ps m ACOM (t A ) \t k 



oc fc-err+r . (231) 



From Fig. [7] we find that the matter power spectrum in the f(R) model (|47|) with n = 3/4 is in fact larger than that 
in the ACDM model on smaller scales. 

The galaxy matter power spectrum is modified by this effect. Meanwhile the CMB spectrum is hardly affected 
except for very large scales (for the multipoles £ = 0(1)) at which the Integrated Sachs- Wolfe (ISW) effect becomes 
important. Hence there is a difference for the spectral indices of two power spectra, i.e. 



A»(U) = f^. (232) 

For larger n the redshift Zk can be as close as Z\, which means that the estimation (|232[) is not necessarily valid in 
such cases. Moreover the estimation (1232|) does not take into account the evolution of 8 m after z = z\ to the present 
epoch (z = 0). It was found in Ref. )71] that the estimation (|232p agrees well with the numerically obtained An(tji) 
for n < 2. 

In order to discuss the growth rate of matter perturbations, it is customary to introduce the growth index [34 ll ] 

Ho m 

where Ct^ = n 2 p m /(3H 2 ). In the ACDM model 7 is nearly constant for the redshifts < z < 1, i.e. 7 ~ 0.55 
[342l |343| . In f(R) gravity, if the perturbations are in the GR regime (M 2 » k 2 /a 2 ) today, 7 is close to 0.55. On 
the other hand, if the transiti on to the scalar-tensor regime occurs at the redshift 2% larger than 1, the growth index 
tends to be smaller than 0.55 (344L l345l |. Since < fl m < 1, the smaller 7 implies a larger growth rate. 
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For the wave numbers relevant to the linear regime of the matter power spectra (0.01 /iMpc -1 < k < 0.2 /iMpc -1 , 
where h ~ 0.7) the viable f(R) models (T4"5)) - (l4"5)) can give rise to the growth index 70 = 0.40 today. Depending on the 
wave numbers k, 7n can be dispersed in the regime 0.40 < 70 < 0.55, or 70 can show the convergence in the regime 
0-40 < 70 < 0.43 [344, 345|. Moreover, when 70 is small, the growth index exhibits a large variation even for low 
redshifts (0 < z < 1). Although the present observational constraints on 7 are quite weak, the unusual evolution of 7 
can be useful to distinguish the f(R) models from the ACDM in future observations. 

After the system enters the epoch of cosmic acceleration, the wave number k can again become smaller than aM. 
Hence the fc-dependence is not necessarily negligible even for z < z\. However numerical simulations show that An (to) 
is not much different from An(t\) derived by Eq. fl232[) [7lj . Thus the analytic estimation (j232[) is certainly reliable 
to place constraints on model parameters except for n»l. Observationally we do not find any strong difference for 
the slopes of the spectra of LSS and CMB. If we take the mild bound An(t\) < 0.05, we obtain the constraint n > 2 
(69j . In this case the local gravity constraint (|5T1) is also satisfied. 

For th e wave nu mbers k > 0.2 ftMpc -1 we need to take into account the non-linear effect of density per turbation s. 
In Refs. [346l4349j the aut hors carried out A-body simulations for the f(R) model (|46p (see also Refs. [35CH . |351| ). 
Hu and Sawicki (HS) [352] proposed a fitting formula to describe the non-linear power spectrum based on the halo 
model. Koyama et al. [3531 ] studied the validity of the HS fitting formula by comparing it with the results of A-body 
simulations and showed that the HS fitting formula can reproduce the power spectrum in A-body simulations for the 
scales k < 0.5/iMpc -1 . In the quasi non-line ar re gime a normalized skewness, S3 = (<5^) / ' {8 2 n ) 2 , has been evaluated 
in f(R) gravity and in Brans-Dicke theories |354| . The skewness in f(R) dark energy models differs only by a few 
percent relative to the value S3 = 34/7 in the ACDM model. 

The modified growth of matter perturbations also affects the evolution of the gravitational potentials \t and <E>. 
The effective pote ntial if> = $ — \P is important to discuss the ISW effect on the CMB as well as the weak lensing 
observations [355| . For the modes deep inside the Hubble radius, Eq. (|219[) gives 

tl> = 3nm5 m (jp\ . (234) 

Note that for the large-scale modes relevant to the ISW effect in CMB we need to solve the full perturbation equations 
without using the quasi-static approximation on sub-horizon scales. For viable f(R) models, however, numerical 
simulations show tha t the result (|234l) can be trustable even for the wave numbers close to the Hubble radius today 
(i.e. k > a H ) [l55| . In the ACDM model the potential i/j remains constant during the standard matter era, but 
it decays after the system enters the accelerated epoch, producing the ISW contribution for low multipoles on the 
CMB power spectrum. In f(R) gravity the additional growth of matter perturbations in the region z < Zk changes 
the evolution of ip. 

From the observation of the CMB angular power spectrum, the constraint on the deviati on p arameter B = 
(RH/RH)m from the ACDM is weak. The value of B today is constrained to be Bq < 4.3 [356| . There is an- 
other observational constraint coming from the angular correlation between the CMB temperature field and the 
galaxy number density field induced by the ISW effect. The avoidance of a large anti-correlation between the obser- 
vational data of CMB and LSS places an upper bound of Bo < 1. Since this roughy corresponds to m(z = Q) < 1, the 
CMB observations do not provide tight constraints on f(R) models relative to the matter power spectrum of LSS. In 
weak lensing observations, the mo dified evolutio n of the lensing potential ip directly leads to the change even for the 
small-scale shear power spectrum |340ll355ll357l | . Hence this can be a powerful tool to constrain f(R) gravity models 
from future observations. 



B. Brans-Dicke theory 

Let us proceed to discuss the evolution of matter perturbations in Bran-Dicke theory described by the action (|76|) 
with the potential U((t>) and the coupling F(<f>) = e~ 2( 2*. We are mainly interested in large coupling models with \Q\ 
of the order of unity [76|, |8l|, l82| , as this gives rise to significant deviation from the ACDM. We define the field mass 
squared to be M 2 = U^. If the scalar field is light such that the condition M < H is always satisfied irrespective of 
high or low density regions, the coupling Q is constrained to be |Q| < 10 -3 from local gravity constraints. Meanwhile, 
if the mass M in the region of high density is much larger than that on cosmological scales, the model can satisfy local 
gravity constraints under the chameleon mechanism even if |Q| is of the order of unity. Cosmologically the mass M 
can decrease from the past to the present, which can allow the transition from the "GR regime" to the "scalar-tensor 
regime" as it happens in f(R) gravity. An example of the field potential showing this behavior is given by Eq. (I98)) . 

As in f(R) gravity, the matter perturbation S m obeys Eq. (|213[) . The difference appears in the expression of the 
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gravitational potential tf>. In Fourier space the scalar metric perturbations obey the following equations [76l 13251 . |32€ 



p . 1 
— <f> + 3H(HV- $) = 



a 



2F 



uj>6j) + liuj,^ 2 - FjR + 2V t4> ) 



k 2 



3H + 3H Z - — 5F - 3HSF + (3HF - w<^)* + 3F(i/tf - $) + 5 Pri 
a 2 



(235) 



8<j>- 



a 2 
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where <5(/> is the perturbation of the field <p, u = (1 — 6Q 2 )F, and 



2w 



-F.^fl, 



* + $ 



5R = 2 



3($ - i? - 12H(Hty - $) + fc 2 /a 2 - 3i? * + 2(/c7a 2 )$ 



(236) 
(237) 

(238) 



Provided that the field is sufficiently heavy to satisfy the conditions M 2 ^> R, one can employ the approximation 
[(2E/,0 — F ! d,R)/2v],(j> — M 2 /ui in Eq. (|236p . The solution to Eq. (|236l) consists of the sum of the matter-induced mode 
<5</>ind sourced by the matter perturbation and the oscillating mode <50 O sc, i-e. 6<f> — 54>\ n d + S<fi OS c (as in the case of 
f(R) gravity). 

In order to know the evolution of the matter-induced mode we employ the quasi-static approximation on sub-horizon 
scales. Under this approximation, we have 5R ind ~ 2{k 2 /a 2 )[$> - (F <t> /F)S<p- md ] from Eqs. ((235]) and (j238)) . where the 
subscript "ind" represents the matter-induced mode. Then from Eq. (|236[) we find 



find 



2QF 



A: 2 



(fc 2 /a 2 )(l-2Q 2 )F + M 2 a : 



Using Eqs. (g55J) and (j2"3"7) we obtain 



^ m (fc 2 /a 2 )(l + 2Q 2 )F + Af 2 



2F 



{k 2 /a 2 )F + M 2 



^ m (fc 2 /a 2 )(l-2Q 2 )F + Af 2 



2F 



{k 2 /a 2 )F + M 



f2 



(239) 



(240) 



In the massive limit M 2 /F 3> k 2 /a 2 , we recover the standard result of General Relativity. In the massless limit 
M 2 /F < fc 2 /a 2 , it follows that (k 2 /a 2 )^ ~ -(5p m /2F)(l + 2Q 2 ) and (fc 2 /a 2 )$ ~ (fy m /2F)(l - 2Q 2 ). 
Plugging Eq. (|240[) into Eq. (|218[) , we obtain the equation for matter perturbations [76| 



<5 m + 2HS m - 4irG c ffp m 5 m = , 
where the effective gravitational coupling is 

_ G (fc 2 /a 2 )(l + 2Q 2 )F + M 2 
cff ~F (k 2 /a 2 )F + M 2 

We have recovered the bare gravitational constant G. In the massless limit (M 2 <C k 2 /a 2 ) this reduces to 

G . G4 + 2wbd 



G cff ~ -(1 + 2Q 2 



F 3 + 2w BD 



(241) 



(242) 



(243) 



where in the last line we have used the relation (|77[) between the coupling Q and the Brans-Dicke parameter wbd • In 
f(R) gravity we have wbd = and hence G c g = 4G/(3F). 

Let us consider the evolution of the oscillating mode of perturbations. Using Eqs. (|235|) and (|236p for sub-horizon 
modes (k 2 /a 2 » H 2 ), the gravitational potentials can be expressed by S(j) osc (note that 5p m = for the oscillating 
mode). From Eq. (|238[) the perturbation of R corresponding to the oscillating mode is given by 



6R n 



6Q 



+ 3HS<j) osc + (k 2 /a 2 )60 o 



(244) 
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Substituting Eq. (|244l) into Eq. ([236]) . it follows that 

54> osc + 3H6<p osc + (k 2 /a 2 + M 2 /F) 8dp osc ~ , (245) 

which is valid in the regime M 2 R. 

When \Q\ = 0(1) the field potential U{4>) is required to be heavy in the region of high density for the consistency 

with local gravity constraints. We take the potential ([55]) as an example of a viable model. During the matter era the 

i 

field 4> settles down at the instantaneous minima characterized by the condition (llOOp . Then we have that <f> oc pm 1 

2-p 

and M 2 oc pm p during the matter-dominated epoch. The field <j) can initially be heavy to satisfy the condition 
M 2 j F ^> k 2 /a 2 for the modes relevant to the galaxy power spectrum. Depending upon the model parameters and 
the mode k, the mass squared M 2 can be smaller than k 2 /a 2 during the matter era [76j. 

In the regime M 2 /F ^> k 2 /a 2 the matter perturbation equation (|241j) reduces to the standard one in Einstein grav- 
ity, which gives the evolution S m oc t 2 ^ 3 . For the model (|98| the matter- induced mode of the field perturbation evolves 

2(4-p) p / _ i \ 

as 54>i n d oc dp m /M oc t 3 ^~ p > . Meanwhile, the WKB solution to Eq. (|245|) is given by 5(j) osc oc ft 1 -*) cos ( ct I, 

2 

where c is a constant. Since the background field <fi during the matter era evolves as (f) oc t 1 -p , we find 

(5</>/0 = (50 ind + 5<f> osc )/ci) ~ cit 2/3 + c 2 r 3(H) cos (ct - ^) . (246) 

As long as the oscillating mode is initially suppressed relative to the matter-induced mode, the latter remains the 
dominant contribution in the subsequent cosmic expansion history. 

In the regime M 2 /F -C k 2 /a 2 the effective gravitational coupling is given by Eq. (|243p . Solving Eq. (|241[) in this 
case, we obtain the solution 



S m oc t (V»+«g?-i)/6 _ (247) 

Setting Q = — 1/V6, this recovers the solution 6 m oc ^(v^- 1 )/ 6 i n metric f(R) gravity. 

The potential (|98[) has a heavy mass M much larger than H in the deep matter-dominated epoch, but it gradually 
decreases with time. Depending on the modes k, the system crosses the point M 2 /F — k 2 /a 2 at t = % during 

2-p 

the matter-dominated epoch. Since the field mass evolves as M oc t l -p during the matter era, the time tk has a 

_ 3(1— p) 

scale-dependence given by tk oc k 4 -p . Since < p < 1 the smaller scale modes (i.e. larger k) cross the transition 
point earlier. During the matter era the mass squared is approximately given by 

2-p 

f 2 _ , 1 - P ^2 I Pn 



M '* W¥cW^ Q \% 1 (248) 



Using the relation p m — 3F(jQ^ Hq (1 + z) 3 , the critical redshift zu at time tk can be estimated as 

- 1 , (249) 
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a H Qj (1 - p)i-P (3F n^P H 2 



where ao is the scale factor today. The critical redshift increases for larger k/(aoHo) and for smaller p. If k/(aoHo) = 
600 and p = 0.7, Eq. (|2"4^1) gives z k = 3.9. 

Defining the growth rate of matter perturbations as in Eq. (|233[) . it follows that the asymptotic values of f$ in the 
regions t <C tk and t tk are given by fg = 1 and fg — (y/25 + 48Q 2 — l)/4, respectively. Numerical simulations 
show that the growth rate reaches a maximum value around the end of the matter era and then it starts to dec rease 
during the epoch of cosmic acceleration [76] . The observational constraint on fg reported by M cDon ald et al. |358j | 
is fg — 1.46 ± 0.49 around the redshift z = 3, whereas the data reported by Viel and Haehnelt |359] in the redshift 
range 2 < z < 4 show that even the value fg — 2 can be allowed in some of the observations. If we use the criterion 
fg < 2 with the analytic estimation fg = (a/25 + 48Q 2 — l)/4, we obtain the bound Q < 1.08. Note that the growth 
index today can be smaller than 0.4 for \Q\ larger than 0.4, so this will be also useful to place tight bounds on Q in 
future observations. 

The relative difference of the matter power spectrum Pg m at time t — t\ (at which a = 0) from that in the ACDM 
is given by 

„( \Z25 + 48Q a _ 1 2 \ 

Pg m {t A ) ft\\ \ 6 3 J (1-P)(V25 + 48Q 



p£f BM (t A ) U 



oc k ~ . (250) 



40 



In Fig. [7] we plot the matter power spectrum for Q = 0.7 and p = 0.6, which deviates from that in the ACDM model 
on small scales. The estimation (|250[) shows fairly good agreement with numerical results [340j |. 

From Eqs. (|240|) we find that the effective gravitational potential ip = $ — \P satisfies the same equation as ()234|) . 
Since the ISW effect induced by the modified evolution of tp is limited on large-scale CMB perturbations irrelevant 
to the galaxy power spectrum, there is a difference between the spectral indices of the matter power spectrum and of 
the CMB spectrum on the scales, k > 0.01ft Mpc -1 : 

(l-p)(^5 + 48^- 5)j 
4 — p 

This reproduces the result (|232j) in f(R) gravity by setting Q = —l/y/6 and p = 2n + 1. If we use the criterion 
An(t\) < 0.05, as in the case of the f(R) gravity, we obtain the bounds p > 0.957 for Q — 1 and p > 0.855 for 
Q = 0.5. As long as p is close to 1, it is possible to satisfy both cosmological and local gravity constraints for \Q\ < 1. 



C. DGP model 



In this section we study the evolution of linear matter perturbations in the DGP braneworld model. The discussion 
below is valid for the wavelengths larger than the Vainshtein radius r* . For the radius r smaller than r* the non-linear 
effect coming from the brane-bending mode becomes crucially important. The pert urbe d metric in the 5-dimensional 
longitudinal gauge with four scalar metric perturbations "J, $, B, E is given by [87], l36l| 

ds 2 = -(1 + 2#)n(i, yfdt 2 + (1 + 2$)A(t, y) 2 S l:j dx l dx 3 + 2r c B >i dx i dy + (1 + 2E)dy 2 , (252) 

where the brane is located at y = in the 5-th dimension characterized by the coordinate y (we are considering a flat 
FLRW spacetime on the brane). Note that B can be identified as a brane bending mode describing a perturbation 
of the brane location and th at r 9 is the crossover scale defined in Eq. (|130l) . The background solution describing the 
self-accelerating Universe is [218| 



n(t,y) = l + H(l + H/H 2 )y, A(t,y) = o(t)(l + Hy) . (253) 

The Hubble parameter H = a /a satisfies Eq. (1133j) with e = +1. 

In what follows we neglect the terms suppressed by the factor aH/k <C 1 because we are considering sub-horizon 
perturbations. We also ignore the terms such as (A' /A)<f>' , where a prime represents a derivative with respect to y. 
This comes from the fact that $' is of the order of (fc/a)$, as we will show later. The time-derivative terms can be 
also dropped under the quasi-static approximation on sub-horizon scales. Then the perturbed 5-dimensional Einstein 
tensors 5Gg obey the following equations locally in the bulk [36 lj : 

SG° = 3$" + Av 2 * + ^(E- r c B>) - 2^ (*f) V 2 B = , (254) 

5G) : 



A 2 A 2 c ' A 2 \A 



■^(V'Vj - <5*V 2 )($ + * + E - r c B') + <S*(¥" + 2$") + ^(V J V, - ^V 2 ) ( — + ^ J B = , (255) 



SGf = -(*' + 2$'),i = 0, (256) 
SGI = ^ + 2$)-j( 2 | + ^V 2 fi = 0. (257) 

Taking the divergence of the traceless part of Eq. (|255[) , we obtain 

J (ff + * + E r c B>) - ^ + J) V 2 B = . (258) 

The consistency between Eqs. (I256P and (|257l) requires that 

B' = 0, + 2$' = . (259) 

From Eqs. (|257| and (|258|) we find 

V 2 1 V 2 r n' 
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Substituting Eqs (1257) and (|260l) into Eq. ([2541 together with the use of Eq. (|259| . we obtain 

*'' + $*-^V 2 i? = 0. (261) 

Under the sub-horizon approximation (k/aH 3> 1) the solution to Eq. (|261l) . upon the Fourier transformation, is 
given by 



n 

* r r B = 



Cl {l + Hy)- k ' aH + c 2 (l + Hy) k / aH ] , (262) 



where c\ and Ci are integration constants. In order to avoid the divergence of the perturbation in the limit y — > oo 
we choose c 2 = 0. 

The junction condition at th e br ane can be expressed in terms of an extrinsic curvature and an energy- 
momentum tensor on the brane [362] : 

K^ v - Kg^ = -k 2 5) + r c G^ v , (263) 

where K = Kj^. The extrinsic curvature is defined as = h^V\n v , where n v is the unit vector normal to the 
brane and = — n^n^ is the induced metric on the brane. The (0, 0) and spatial components of the junction 
condition (I263p give 



-V 2 $ = -K 2 {i) 6 Pm + —V 2 B - , (264) 



$ + # = (265) 
+ 2$' = , (266) 



where Sp m is the matter perturbation on the brane. Equation (|266[) is consistent with the latter of Eq. (|259|) . 

From Eq. (I262[) it follows that $' ~ (k/a)Q in Fourier space. For the perturbations whose wavelengths are much 
smaller than the cross-over scale r c , i.e., r c k/a ^S> 1, the term (3/r c )<&' in Eq. (|264[) is much smaller than (fc 2 /a 2 )$. 
In Fourier space Eq. (|264|) is approximately given by 

2p k 2 

—<S> = K 2 {A) 5p m + ^B. (267) 

Using the projection of Eq. f|257[) as well as Eqs. (|265[) and (|267l) . we find that metric perturbations * and $ obey 
the following equations 

a 2 

where 



^* = 1 1 + yj) s '- ?* = ^l 1 -3?J ,A "' (268) 



«'> si -| £ (4 + ^)= i - 2 ^( i+ 3^)- « 

The matter perturbation <5 m satisfies the same form of equation as given in (|218[) for the modes deep inside the horizon 
[360j |. Substi tutin g the former of Eq. (|268l) into Eq. (|2 18[) . we find that the matter perturbation obeys the following 
equation [H HI 

5 m + 2HS m - ATrG cS p m S m = , (270) 

where 

G«=(l + ±)o. (271) 

From Eq. f|268[) the effective gravitational potential ip = $ — obeys the same equation as (|234[) . 

In the deep matter era one has Hr c 3> 1 and hence f3 ~ —Hr c , so that /3 is largely negative (\f3\ ^> 1). In this regime 
the evolution of the matter perturbation is similar to that in General Relativity (S m oc i 2 / 3 ). The solutions finally 
approach the de Sitter attractor characterized by H^s = l/r c . At the e Sitter solution one has {3 ~ 1 — 2Hr c ~ — 1. 
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Since 1 + 1/ (3/3) ~ 2/3 , the growth rate in this regime is smaller than that in GR. The growth index is approximately 
given by 7 w 0.68 343], which is different from the value 7 ~ 0.55 in the ACDM model. If the future imaging survey 
of ga laxies can constrain 7 within 20 %, it may be possible to distinguish the DGP model from the ACDM model 
(363j . 

Comparing Eq. (|271l) with the effective gravitational coupling (|243[) in Brans-Dicke theory with a massless limit (or 
the absence of the field potential), we find that the Brans-Dicke parameter wbd has the following relation with /3: 

wbd = ~ (/3- 1). (272) 

Since j3 < for the self-accelerating DGP solution, this implies that cjbd < —3/2. Since in this case the kinetic energy 
of a scalar field degree of freedom is negative in the Einstein frame, the DGP model contains a ghost mode. The 
solution in another branch of the DGP model is not plagued by this problem, because the minus sign of Eq. (|269[) is 
replaced by the plus sign. The self accelerating solution in the original DGP model can be realized at the expense of 
an appearance of the ghost state. 



D. Galileon gravity 

1. Covariant Galileon gravity 

In co varia nt Galileon gravity described by the action (|140|) the evolution of matter density perturbations was studied 
in Ref. [364{ . In spite of the complexities of full perturbation equations, they are simplified under the quasi-static 
approximation on sub-horizon scales. Under this approximation the matter perturbation obeys the same equation as 



(12701) with a different effective gravitational coupling G c g. In Ref. 364] it was found that G B s is independent of the 
wave number k, as in the DGP model. G e g is close to the gravitational constant G in the asymptotic past, but it 
starts to deviate at the late cosmological epoch. 

We define the effective gravitational potential ^ = $ — as well as the anisotropic parameter 77 = — $/W. Under 
the quasi-static approximation on sub-horizon scales we obtain [364j 

^3%^^X) 2 ' (273) 

where f2 m and 5 m and the density parameter and the perturbation of non-rclativistic matter, respectively. In f(R) 
gravity, Brans-Dicke theory, and the DGP model the effective gravitational potential obeys Eq. (|234[) for the modes 
deep inside the Hubble radius. In Galileon gravity the combination (G e ff/G)(l +rf)/2 is different from 1. This means 
that the effective gravitational potential may acquire some additional growth compared to other models. In Galileon 
cosmology there are three different regimes characterized by (i) r\ <c 1, 7*2 <C 1, (ii) T\ = 1, Ti *C 1, and (i ii) n = 1, 
r- 2 = 1, where r\ and r 2 are defined in Eq. (|150|) . In these regimes we can estimate G e ff and 77 as follows 364]. We 
stress that these analytic results are valid for the modes deep inside the Hubble radius. 

• (i) n < 1, r 2 < 1 

Expanding G e ff and 77 about n = 0,r 2 — 0, it follows that 

G off 1 , {255 n ,211 \ , ml 2 f 12g a 589 \ 2 



G l+^—/3 + —a ri jr 2 + 0(r A 2 ), V = 1 + I — + —an j r 2 + 0(r 2 2 ) , (274) 

where a and /3 are defined by Eq. (|149[) . Since [3 > to avoid ghosts (for the branch r 2 > 0), we have G e g > G 
and 77 > 1 in this regime. This means that the growth rates of S m and ip are larger than those in the ACDM 
model. 

• (ii) n = 1, T2 <C 1 
Expansion of G e ff and 77 about r 2 = gives 

Geff . , 291a 2 + 702/3 2 - 933a/3 + 20a - 84/3 + 4 a 

~G = 1 + 2(10a-9/3 + 8) r2 + ' (275) 

, 3(126a 2 + 306/3 2 - 405a/3 + 4a - 30/3) m/ 2 , 
" = 1 - J 2(10a-9/3 + 8) ^ 2 + ° (r2) ' (276) 

The evolution of G c ff and 77 depends on both a and j3. If a = 1.4 and (3 = 0.4, for example, we have 
Geff/G ~ 1 + 4.31r2 and 77 ~ 1 — 5.11r2, respectively. In this case G c g > G, but 77 is smaller than 1. 
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Figure 8: Evolution of the perturbations for a — 1.37 and /3 = 0.44 with the background initial conditions n = 0.03 and 
T2 = 0.003. (Left) 8 m /a versus z for the wave numbers (a) k — 300aoi/o, (b) k — 30aoHo, and (c) k = 5aoHo- (Right) ip versus 
z for the wave numbers (a) k = 300ao-Ho, (b) k = 10ao-Ho, and (c) k = 5aoHo- Note that 8 m /a and ip are divided by their 
initial amplitudes 5 m {ti) / a(ti) and ip{ti), respectively, so that their initial values are normalized t o be 1. The bold dotted lines 
show the results obtained under the quasi-static approximation on sub-horizon scales. From Ref. |364l ]. 



• (iii) n = 1, r 2 = 1 
At the dS point we have 



G e ff 1 



G 3(a-2/3) 
which means that there is no anisotropic stress. 



(277) 



Recall that the late-time tracking solutions are favored from observational constraints at the background level. In 
this case the solutions start from the regime (i) and finally approach the de Sitter fixed point with a short period of 
the regime (ii). In Fig. |H] we plot S m /a and ip versus the redshift z for a = 1.37 and f3 — 0.44 with the background 
initial conditions ri = 0.03 and r 2 = 0.003. In this case the solutions approach the tracker at late times. The initial 
conditions of perturbations are chosen to recover the GR behavior in the asymptotic past. Note that these results 
are obtained by numerical integration of the full perturbation equations without quasi-static approximations. For 
the mode k — 300ao-ffo the numerical result shows excellent agreement with that obtained under the quasi-static 
approximation on sub-horizon scales. The difference starts to appear for the modes fc/(ao-ffo) < O(10). From the left 
panel of Fig. [8] we find that, on larger scales, the growth of 8 m tends to be less significant. For the modes k a^Ho 
the matter perturbation evolves faster than a during the matter era. 

From the right panel of Fig. [S]we find that, unlike the ACDM model, ip changes in time even during the matter era 
for the modes k 3> aH. Before reaching the tracker we have G c s/G ~ 1 + 255/3^/8 > 1 and 77 ~ 1 + 129/3^/8 > 1 from 
Eq. (|274l) . Then the growth rates of ip and S m get larger than those in GR. In particular the term (G c g/G)(l +rj)/2 
in Eq. (|273[) is larger than 1, which leads to the additional growth of ip to that coming from S m . In Galileon gravity 
the unusual behavior of the anisotropic parameter rj leads to the non-trivial evolution of perturbations. For the model 
parameters a — 1.37 and (3 — 0.44, Eq. (|277|) gives G c s ~ 0.68G at the de Sitter fixed point. Since in this case G e s 
is smaller than G, ip begins to decrease at some point after the matter era. 

For the large-scale modes relevant to the ISW effect in CMB anisotropics (k/(aoHo) < 10), ip is nearly constant in 
the early matter-dominated epoch. However, as we see in Fig. [51 ip exhibits temporal growth during the transition 
from the matter era to the epoch of cosmic acceleration. The characteristic variation of ip in the Galileon model may 
leave interesting observational signatures on the large-scale CMB anisotropies. 

For the model pa rame ters constrained by SN la (Union2)+CMB+BAO data sets, i.e. a = 1.404 ± 0.057 and 
/3 = 0.419 ± 0.023 [259| | . the effective gravitational coupling at the de Sitter solution is restricted in the range 
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0.5G < G G ff < 0.72G. If the model parameters are close to the upper limit a — 2(3 + 2/3 of the allowed parameter 
space at the background level (i.e. G D g is close to 0.5G at the de Sitter point), the parameter r\ tends to show a 
divergence during the transition from the matter era to the epoch of cosmic acceleration. If G c g is larger than 0.66G, 
we find that such divergent behavior is typically avoided. For the viable model parameters the evolution of 8 m and ip 
exhibits more or less the similar property to that shown in Fig. [8] 



2. Modified Galilean gravity 



Finally we study modified Galileon theories in which the term D<p(d ^(pd^ 1 <p) is generalized to ^((p)U\(p(d ^<pd^ (p) . Let 



us consider general theories described by the action 



S = d 4 X\/— g 



d 4 x£ 



M ■ 



where 



f{R,^X) = h{R,cj>)+f 2 {4>,X) 



(278) 



(279) 



We introduce two mass scales associates with the field <p and the scalar gravitational degree of freedom, respectively, 
as 



Mi 



-SmI* , M% = F/(3F R ) , 



(280) 



where F= df/dR. 

The full perturbation equations for the perturbed metric (|207[) are given in Ref. [183]. Under the quasi-static ap- 
proximation on sub-horizon scales the matter perturbation satisfies Eq. (|270[) with the effective gravitational coupling 

|i~8l 



G e ff = 



1 1+481 

8irF 1 + 3si 



[J> + 2(l + 4 Sl )^ 2 
(1 + 4si)/xF 



(281) 



where 



H = (1 + 3ai)(f, x + 2s 2 ) + 3F%/F + 2£[4(1 + 3«i)$ + 2Hj>) - 2F^ 2 /F - 2^ 4 (1 + 4s x )/F] , (282) 

2//•Q„2n / r2^ „ — „2n / /-2/;„2 (283) 



s x = k 2 /{3a 2 M 2 R ) , s 2 = a 2 Ml/k 2 . 
The anisotropic parameter 77 = — $/\E r is given by 

(1 + 2«i)(/, x + 2s 2 ) + 2^/F + 4£[2(1 + 2s x )(0 + 2ff0) - F^ip 2 jF 



V 



(1 + 4ai)(/ iJr + 2s 2 ) + 4F 2 /F + 8£(1 + 4«i)(0 + 2#» 
The effective gravitational potential ip = $ — \& obeys the equation of the form (I273|) . whose explicit form is 



(284) 



ip = 3Q m 5 m — — 



a_H\ 2 (1 + 3si)[/, x + 2s 2 + 8!j(4> + 2H<j>)] + 3F|/F - 2£F^ 2 /F 

k ) (1 + 3 Sl )[J> + 2s 2 + 8£(0 + 2ff$] + 3F 2 ,/F - ^[F^/F + £</> 4 (l + 4 3l )/F] 



(285) 



where fi m = p m /(3FH 2 ). 
In the massive limits M 2 



and M| 



00, i.e. s\ — > and s 2 —> 00 we recover the standard results in 



GR: G ofI ~ 1/(8ttF), r, ~ 1, and V ^ 30 m (5 m (aH/k). The difference appears in the regimes characterized by the 
conditions si > 1 and s 2 < 1, as it happens for the late-time evolution of perturbations in f(R) gravity and in 
Brans-Dicke theory. 

For the theories with £ = 0, Eq. (|285[) gives the standard relation ip — 3Q m 8 m (aH/k) 2 . If £ ^ 0, then this relation 
no longer holds. In this case ip is not directly related with S m due to the additional contribution from the ^-dependent 
term. This can be regarded as the main reas on fo r the anti-correlation between the ISW effect in CMB and the 
large-scale stru cture found for the model (I182[) [239] . 

In Ref. jl83| the authors derived conditions for avoiding the appearance of ghosts and Laplacian instabilities. 
Provided F > the tensor ghosts do not appear. Since c 2 ^ = 1 for the theories (|278p . tensor perturbations have 
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no Laplacian instabilities. For the theories with f rr 7^ the conditions for the avoidance of ghosts and Laplacian 
instabilities of scalar perturbations are given, respectively, by 

24£ii> - 8£,^ 2 + /, x + f,xx4> 2 - F 2 x <j> 2 /F R > , (286) 
f.x + 8(0 + 2Hj>)ti - 16£ 2 </> 4 /(3F) > . (287) 

Similar conditions have been also derived for the theories with f rr ^ 0. The dark energy models based on the action 
(|278p need to be constructed to satisfy these conditions. 



E. Observables in weak lensing 



We have shown that modified gravity models of dark energy generally lead to changes for the growth rate of matter 
perturbations compared to the ACDM model. Since there are two free functions that determine the first-order metrics 
^ and dark energy models can be classified according to how the gravitational potentials are linked to 8 m . In order 
to quantify this, we introduce two quantities q(k,t) and C(k,t) defined by 

{k 2 /a 2 )$ = 4TrGq5 mPm , (288) 
($ + *)/$ = C, (289) 

where G is the 4-dimensional bare gravitational constant. The ACDM model corresponds to q = 1 and £ = (note 
that the cosmological constant does not cluster). Modified gravity models give rise to different values of q and £ relative 
to those in the ACDM model. Therefore the functions q and ( characterize gravitational theories for first-order scalar 
perturbations on small scales. 

In Brans-Dicke theory discussed in Sec. I VII B[ the gravitational potentials are given by Eq. (|240|) on sub-horizon 
scales. In this case we have 

_ 1 (fc 2 /a 2 )(l-2Q 2 )F + M 2 _ 4F(fc 2 /a 2 )Q 2 

q F (k 2 /a 2 )F + M 2 ' ^ (fc 2 /a 2 )(l - 2Q 2 )F + M 2 ' 1 ' 

where we have used the unit 8irG — 1. In the regime M 2 / F 3> k 2 /a 2 (and F ~ 1) it follows that q ~ 1 and £ ~ 0. 
In the regime M 2 / F <C k 2 /a 2 we have q ~ (1 — 2Q 2 )/F and ( ~ — 4Q 2 /(1 — 2Q 2 ), so that the deviation from the 
ACDM model appears. The expression (|290l) covers the case of f(R) gravity by setting Q = — 1/V6. 
In the DGP model the gravitational potentials obey Eq. (|268p and hence 

9=1-1/(3/3), C= 2/(1 -3/3). (291) 

In the deep matter era one has \[3\ ^> 1, so that q ~ 1 and ( ~ 0. The deviation from (q, Q — (1, 0) appears when \(3\ 
decreases to the order of unity, i.e., when the Universe enters the epoch of cosmic acceleration. 

In order to confront dark energy models with the observations of weak lensing, it may be convenient to introduce 
the following quantity [355j 

E = g(l-C/2). (292) 

From Eqs. (|288|) and (I289j) we find that the effective gravitational potential i/j = $ — $ associated with weak lensing 
observations can be expressed as 

if) = 87rG , (a 2 //c 2 )p„ l (5„ l I] . (293) 

In the DGP model and in Brans-Dicke theory we have E = 1 and E = 1/F, respectively. In (modified) Galileon 
theories the term E is of more complicated forms, see e.g., Eq. ([285]) . 

The effect of modified gravity theories manifests itself in weak lensing observations in at least two ways. One is the 
multiplication of the term E on the r.h.s. of Eq. (|293[) . Another is the modification of the evolution of S m . The latter 
depends on two parameters q and £, or equivalently, S and £. Thus two parameters (E, £) will be useful to detect 
signatures of modified gravity theories from future surveys of weak lensing. 
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VIII. CONCLUSIONS 

We have reviewed modified gravitational models of dark energy responsible for the cosmic acceleration today. In 
addition to cosmological constraints such as the presence of a matter era followed by a stable de Sitter solution, we 
require that the models satisfy local gravity constraints. There are two mechanisms for the recovery of GR behavior 
in the regions of high density. The first one is the chameleon mechanism in which the mass of a scalar-held degree 
of freedom depends on the matter density in the surrounding environment. The chameleon mechanism can be at 
work in f(R) gravity and Brans-Dicke theory, as long as the field potential is designed to have a large mass in the 
regions of high density. The second one is the Vainshtein mechanism in which the nonlinear effect of scalar-field self 
interactions leads to the recovery of GR at small distances. This can be applied to the DGP braneworld model and 
Galileon gravity. 

The modified gravity models can give rise to the phantom equation of state of dark energy (wde < — 1) without 
having ghosts, tachyons, and Laplacian instabilities. The deviation of u>de from that in the ACDM model (wde = — 1) 
is not so significant in the models based on f(R) gravity and Brans-Dicke theory, so these models can be compatible 
with the observational constraints at the background cosmology fairly easily. On the other hand, the tracker solution 
in covariant Galileon gravity, which has u>de = ~~ 2 during the matter era, is disfavored from the joint data analysis 
of SN la, BAO, and CMB shift parameters. However the late-time tracking without a significant deviation from 
iode = — 1 is allowed observationally. 

In order to confront the modified gravity models with the observations of large-scale structure, CMB, and weak 
lensing, we have also discussed the evolution of matter density perturbations. In the models based on f(R) and 
Brans-Dicke theories there is a "General Relativistic" regime in which the field is heavy such that M 2 ^> k 2 /a 2 . 
At late times this is followed by a "scalar-tensor" regime (M 2 <C k 2 /a 2 ) in which the gravitational law is modified 
from that in General Relativity. In Brans-Dicke theory the evolution of matter perturbations during the matter era 

changes from S rn oc t 2 / 3 to S m oc t(V 25 + 48< 3 2-1 )/ 6 , where Q is related to the Brans-Dicke parameter via the relation 
3 + 2ojbd = 1/(2Q 2 ). The effective gravitational couplings G e ff in the DGP and Galileon models are independent 
of the wave numbers k. This reflects the fact that the field is massless in those models. In the DGP model the 
growth rate of 6 m is smaller than that in the ACDM, which is associated with the appearance of ghosts. In (modified) 
Galileon gravity, a non-trivial relation between the effective gravitational potential ip and the matter perturbation S m 
leads to the extra growth of ip. 

We summarize the current status of each modified gravity model of dark energy. 

• (i) In f(R) gravity there are some viable dark energy models such as (|45[). (|46|) , and (|48|) , which can be consistent 
with both cosmological and local gravity constraints. For the models (|4"5)l and (|46l) the local gravity constraints 
can be satisfied for n > 0.9 under the chameleon mechanism. If we use the criterion that the difference between 
the spectral indices between the matter power spectrum and the CMB spectrum is smaller than 0.05, then we 
obtain the bound n > 2. In these models the initial conditions of perturbations need to be chosen such that 
the oscillating mode does not dominate over the matter-induced mode in the early Universe. This is associated 
with a weak singularity problem about the divergence of the mass squared M 2 ~ F/(3f i im) for R — » oo, but it 
can be circumvented by including higher-curvature terms such as R 2 to the Lagrangian. 

• (ii) In Brans-Dicke theory, even if the coupling Q is of the order of unity, it is possible to design field potentials 
consistent with both cosmological and local gravity constraints. One of the examples is the potential (|98l) . 
which is motivated by the f(R) models (|4"5|) and (|4"6]) . Depending on the couplings Q the growth of matter 
perturbations is different (Q = — 1/\/6 in metric f{R) gravity). From the observational constraints on the 
growth rate of 8 m , \Q\ is required to be smaller than the order of 1. 

• (iii) In the DGP model the self-acceleration is realized through the gravitational leakage to the extra dimension, 
but the joint analysis using the data of SN la, BAO, and CMB shift parameters shows that the model is in 
tension with observations. Moreover the linear perturbation theory beyond the Vainshtein radius shows that 
the model contains a ghost mode with the effective Brans-Dicke parameter wbd smaller than —3/2. However, 
the modification of the DGP model like Cascading gravity can alleviate this problem. 

• (iv) In covariant Galileon gravity there is a tracker that attracts solutions with different initial conditions to a 
common trajectory. The joint observational constraints at the background cosmology shows that the late-time 
tracking solutions are favored from the data. Cosmological perturbations in (modified) Galileon theories exhibit 
peculiar features because of non-trivial relations between the effective gravitational potential ip and the matter 
perturbation 8 m - 

• (v) The dark energy models described by a Gauss-Bonnet term with a scalar coupling F{<p)Q do not satisfy 
both cosmological and local gravity constraints. The generalized Gauss-Bonnet model in which the Lagrangian 
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density is given by R/2 + f(G) is plagued by a serious problem of the Laplacian instability in the presence of 
matter fluids. There are some viable Lorentz-violating models of dark energy in which the phantom equation 
of state can be realized without having ghosts, tachyons, and Laplacian instabilities. 

We hope to find some signatures for the modification of gravity in future high-precision observations. This will shed 
new light on the nature of dark energy. 
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